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Abstract 


Part  I  -  The  Geometric  Theory  of  Diffraction 

In  Part  I,  a  georaetirlc  method  is  given  for  finding  the  field  in  inhomo- 
geneous  media  containing  smooth  convex  bodies.  The  field  due  to  a  plane  wave  in 
an  unbounded  medium  is  constructed  by  introducing  complex  rays  in  the  refraction 
shadow.  By  extending  Fermat's  principle,  we  explain  the  occurrence  of  certain 
diffracted  rays  in  the  boundary  problems.  By  then  modifying  the  law  of  con- 
servation of  energy,  we  obtain  the  field  along  the  diffracted  rays  in  addition 
to  the  field  in  the  geometric  lit  region.  Certain  diffraction  coefficients 
and  decay  exponents  are  introduced  and  general  formulas  are  obtained  for  them. 
The  theory  is  then  applied  to  several  problems  in  which  the  mediiim  is  plane  or 
cylindrically  stratified  and  the  boundary  planar  or  circular. 

Part  IT  -  Asymptotic  Expansion  of  Solutions  of  Boundary  Value  Problems 

In  Part  II,  we  determine  expressions  for  the  exact  solution  to  various 
boundary  value  problems  corresponding  to  the  special  problems  of  Part  I. 
Asymptotic  forms  are  obtained  by  using  the  method  of  stationary  phase  in  con- 
junction with  the  '/JKB  method  and  Watson  transformations.  In  all  cases,  we 
compare  the  results  df  Part  I  with  these  asjmptotic  expansions  and  find  that 
the  two  are  in  agreement. 
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Part  I,  The  Geometrical  Diffraction  Theory 
1,  Introduction 

In  the  analysis  of  wave  propagation  in  inhomogeneous  media  the  method 
of  geometrical  optics  is  often  employed.  Not  only  is  it  used  to  determine 
where  waves  propagate,  but  also  to  evaluate  the  field  quantitatively.  However 
it  fails  to  yield  any  field  other  than  zero  in  shadow  regions,  i.e.,  regions 
devoid  of  rays.  It  is  our  purpose  to  show  that  the  geometrical  theory  of 
diffraction'-  -  ,  a  recent  extension  of  geometrical  optics,  is  capable  of 
accurately  describing  the  diffracted  field  in  shadows.  It  also  yields 
diffraction  corrections  to  the  geometrical  optics  field  in  illuminated  regions. 

This  theory  differs  from  ordinary  geometrical  optics  by  the  inclusion  of 
certain  additional  rays  called  diffracted  rajs   and  others  called  complex  or 
imaginary  rays.  The  former  are  associated  with  the  boundaries  or  discontinuities 
of  the  medium  while  the  latter  are  associated  with  refraction  shadows  which  occur 
in  the  absence  of  boundaries.  In  addition  to  characterizing  these  new  rays,  the 
theory  also  provides  a  method  for  calculating  the  field  associated  with  them. 

In  Part  I  of  this  article  we  explain  the  theory  and  apply  it  to  a  number 
of  problems.  Then  in  Part  II  we  reconsider  some  of  these  same  problems  from 
the  viewpoint  of  botmdary  value  problems  for  partial  differential  equations.  We 
choose  cases  which  can  be  solved  exactly  by  the  method  of  separation  of  variables. 
Then  we  evaluate  these  solutions  asymptotically  for  large  values  of  ka 
(k  =  -«—  where  X  is  the  wavelength,  a  =  tjrpical  dimension  in  the  problem).  Upon 
comparing  the  results  of  Parts  I  and  II  we  verify  that  the  geometrical  theory 
is  indeed  correct,  at  least  in  these  cases. 

The  geometrical  theory  of  diffraction  has  previously  been  applied  to 
diffraction  by  convex  cylinders^  ■■ ,  diffraction  by  apertures  in  thin  screens '--'•' , 
diffraction  by  smooth  opaque  objects  '-  ■' ,  diffraction  in  waveguides  ^    -' ,  and  other 
problems  L  -■.   In  all  cases  the  results  of  this  theory  coincided  with  the 
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asymptotic  expansions  for  large  ka,  of  exact  solutions  of  boundary  value 
problems,  when  these  were  available  for  comparison.  The  results  were  also 
accurate  for  wavelengths  as  large  as  the  dimensions  of  the  scatterer. 

2»  Ray  tracing 

In  geometrical  optics  a  medium  is  characterized  try  its  index  of  refraction 
n(x).  In  terms  of  n(x)  a  ray  is  defined  as  a  curve  x(t)  which  satisfies  the 
equation 

(1)  *i?  -  V(nV2) 

The  parameter  t  is  normalized  by  the  condition 

(2)  k^  -r? 

2 
It  is  to  be  noted  that  x  denotes  a  vector  so  that  x  denotes  the  square  of  t  he 

length  of  X  .  When  a  ray  hits  an  interface  or  boundary  it  gives  rise  to 

reflected  and  refracted  rays.  The  initial  directions  of  these  rays  are 

determined  by  the  laws  of  reflection  and  refraction,  and  their  subsequent  paths 

are  determined  by  (l)  and  (2). 

The  rays  described  above,  and  the  multiply  reflected  and  refraci;ea  rays 

to  which  they  give  rise,  constitute  all  the  rays  of  geometrical  optics.  In  the 

geometrical  theory  of  diffraction  we  introduce  additional  rays  of  various  kinds 

called  diffracted  rs^s.  The  only  type  of  diffracted  ray  which  we  will  consider 

here  is  the  kind  produced  ^en  a  ray  ia  tangent  to  a  boundary  surface.  We  assume 

that  such  a  ray  splits  at  the  point  of  tangency.  One  part  of  the  ray  continues 

along  the  original  path  of  the  incident  ray.  The  other  part  proceeds  along  the 

boundary  as  a  surface  ray.  At  each  point  the  surface  ray  splits  again, shedding 

a  ray  which  leaves  the  surface  in  the  direction  of  the  tangent  to  the  surface  ray. 

This  shed  ray  then  proceeds  in  accordance  with  (1)  and  (2),  The  path  of  the 
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diffracted  ray  along  the  surface  is  an  optical  geodesic.  By  this  we  mean  that 

the  surface  ray  between  two  points  Q,  and  P,  on  the  surface  is  a  curve  which 

P 
renders  stationary  the  optical  length  I   n[x(s)Jds  among  all  curves  on  the 


(  ^  n[x(s)J< 
^^1 


surface.     Here    s    denotes  arclength. 

The  above  characterization  of  a  diffracted  ray  is  equivalent  to  the 
following  modified  form  of  Ferraat's  principle;     A  diffracted  ray  connecting 


two  points  Q  and  P  is  a  curve  which  renders  stationary  the  optical  length 

P 

n[x(s)]ds  among  all  curves  joining  Q  and  P  and  having  an  arc  on  the  boundary. 

Q 
From  this  principle  it  follows  that  the  diffracted  ray  consists  of  two  segments 


[ 


of  ordinary  rays  which  satisfy   (1)  and  (2),   one  from  "^  tangent  to  the  boundary 
and  another  from  P  tangent  to  the  boundary.     The  two  points  of  tangency  are 
connected  by  a  surface  geodesic  ■v^iich  is  tangent  to  the  two  ordinary  ray  segments. 
These  are  just  the  properties  of  diffracted  rays  described  above. 

In  determining  a  surface  geodesic  the  index  of  refraction  associated  with 
one  of  the  two  sides  of  the  surface  must  be  chosen.     We  will  consider,   for 
simplicity,   only  those  surfaces  which  bound  the  region  within  which  propagation 
occurs.     Therefore  n(x)  will  be  defined  on  only  one  side  of  the  surface,   so  no 
ambiguity  can  arise. 

In  addition  to  the  diffracted  rays,   we  >dll  also  consider  complex  or  imaginary 
rays.     These  are  defined  to  be  complex  valued  solutions  x(t)  of  the  ray  equations 
(1)  and  (2).     Such  rays  can  be  defined  vdien  n(x)  is  an  analytic  function  of  x, 
and  also  in  certain  cases,  \rtien  n  is  piece-wise  analytic.     The  particular  complex 
rays  vrtiich  arise  in  any  problem  are  determined  by  the  incident  rays.     This 
will  become  clear  in  Section  U,  where  these  rays  are  utilized. 

3«     Field  calculation 

In  the  present  theory,   as  in  geanetrical  optics,   a  field  u(P)  is  associated 
with  each  point  P  on  a  ray.     This  field  is  conposed  of  a  phase  (?(s)  and  an 
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amplitnde  A(s),  both  of  vjtiich  may  vary  with  arclength  s  along  the  ray.     The 
phase  is  given  in  terms  of  optical  length  by  the  formula 


(3) 


^{3)  "  9^Q  +  I     n[x(sr|ds 


Here  <j(    is  the  phase  at  the  point  Q  from  which  arclength  s  along  the  ray  is 
measured.     The  direction  of  increasing  s  must  be  the  direction  of  propagation 
along  the  ray. 

The  amplitude  A(s)  is  determined  by  the  principle  of  conservation  of 
energy.     The  optical  form  of  this  principle  states  that  in  a  time  periodic 
(harmonic)  field  the  energy  flux  is  the  same  through  all  cross  sections  of  a 
tube  of  rays.     Upon  applying  this  principle  to  a  narrow  tube  of  rays  containing 
the  ray  from  Q  to  P  we  obtain 


(U) 


A^   (s)  n(s)  dCT"  A^  n     dcT 

000 


area  of  the  tube  at  P,  while  A  ,   n    and  dcr   are  the  corresponding  quantities 

000 


Here  A,   n  and  dcr  denote  the  amplitude,   index  of  refraction  and  cross- sectional 
area  of  the  tube  at  P,  whi 
at  Q.     From   (U)  we  obtain 

(5)  A(s)  -  A 


n  do''1 

o  o 

n(s)     da- 


1/2 


The  ratio  do' /dcr  is  the  limit  of  the  area  ratios  as  the  tube  shrinks  to  the  ray. 
By  combining  (3)  and  (5)  we  find  for  the  field  at  P 


(6) 


u(P)  «  A 


n 
o 


dcr" 
o 

do- 


e:qp  ik^(2f^  +         n[x(sr]dsl 


In  this  equation  k  ■  ca/c  is  the  propagation  constant  defined  in  terms  of  co,  the 
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angular  frequency  of  the  field  and  c,the  reference  propagation  speed  used  in 
defining  n.  A  reference  wavelength  X  is  related  to  k  by  the  equation  k  »  2n/X. 
The  field  (6)  is  associated  with  one  ray  through  P.  The  total  field  at  P  is 
defined  to  be  the  svm   of  the  fields  associated  with  all  the  rays  through  P. 
We  denote  by  u,,  u  ,  u.  and  u,  respectively  the  fields  associated  with  incident 
rays  (i.e.,  rays  from  the  source),  reflected  rays,  transmitted  (refracted) rays 
and  diffracted  rays.  We  also  call  u.  +  u  +  u.  the  geometrical  optics  field  u  . 
Thus  the  total  field  is 

(7)  u(P)  -  u  (P)  +  u  (P)  . 

The  prescription  (6)  for  the  calculaticn  of  the  field  must  be  supplemented 
to  be  used  on  reflected  and  refracted  rays.  On  these  rays  we  assume  that  the 
initial  value  of  the  field  at  the  point  of  reflection  Q  is  proportional  to  the 
incident  field  at  Q.  Thus  u  (Q)  -  Ru.(Q),  u.(Q)  =  Tu.(Q).  The  proportionality 
factors  R  and  T  are  called  the  reflection  and  transmission  coefficients.  We  assume 
that  they  depend  upon  the  nature  of  the  field,  the  local  properties  of  the  media  at 
the  point  of  reflection  and  the  angle  of  incidence.  Because  they  depend  only 
upon  local  properties,  they  can  be  determined  from  the  analysis  of  plane  wave 
reflection  and  transmission  by  a  plane  interface.  Once  they  are  known,  u  can 
be  calculated  completely  for  any  given  u. .  As  an  example,  if  the  reflecting 
surface  is  characterized  by  an  impedance  Z  then 

(8)  R  =  [n  cos  6  -  z]  [n  cos  6  +  z]"  . 

Here  n  and  5  are  the  index  of  refraction  and  the  angle  of  incidence  at  the  point 
of  reflection.  This  result  can  be  obtained  by  inserting  the  expressions  of  the 
fonn  (6)  for  ti.  and  u  into  the  impedance  boundary  condition  Bvi/dv   =  ikZu.  In 
the  limiting  cases  Z  »  0  and  Z  =  oo,  (8)  yields  R  »  1  and  R  =  -1  respectively. 
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Equation  (6)  is  also  applicable  to  diffracted  rays  with  the  exception  of 
those  portions  of  these  rays  which  lie  on  the  surface.  Equation  (3)  for  the 
phase  is  still  applicable  on  these  portions,  but  (5)  must  be  modified.  To 
modify  it  first  denote  by  Q  the  point  of  tangency  of  an  incident  ray  with  the 
surface.  Then  u,(Q,)  is  the  field  on  the  diffracted  (or  surface)  ray  at  Q^.  We 
now  assume  that  the  diffracted  field  at  Q,  is  proportional  to  the  incident  field 
there 

(9)  M^iQ^)   -  D(Q3^)u^(aL). 

The  proportionality  factor  D((i. )  we  call  the  diffraction  coefficient.  We  assvune 
that  it  depends  upon  the  nature  of  the  field,  the  properties  of  the  media  and  the 
shape  of  the  surface  at  Q^,  and  the  propagation  constant  k.  In  the  case  of  vector 
fields  D  is  a  matrix.  Before  describing  how  D  can  be  deteimined  we  will  consider 
the  variation  of  the  diffracted  field  along  the  diffracted  ray. 

Since  the  phase  of  the  diffracted  field  is  given  by  (3),  we  need  consider 
only  the  amplitude.  According  to  (3),  the  phase  of  u,(Q,)  equals  that  of  u.(Q,). 
Therefore  we  have  from  (9) 

(10)  /^(Q^)  =  D(Q^)A^(Q3^). 

Let  us  denote  by  t  the  distance  from  Q.  along  the  diffracted  ray  and  by  -'^^(t) 
the  amplitude  of  the  diffracted  field  at  t.  To  determine  -A^(t)  we  will  consider 
a  narrow  band  of  rays  on  the  surface  and  apply  to  it  the  energy  conservation 
principle.  We  first  denote  by  do-(t)  the  width  of  the  band.  Then  the  energy  fltuc 
along  the  band  is  proportional  to  ^ndo'.  We  assume  that  this  flux  diminishes  due  to 
radiation  via  the  shedding  of  diffracted  rays  from  the  band.  The  energy  removed 
by  the  shed  rays  must  be  proportional  to  the  energy  carried  in  the  band  so  we 
assume  that 


(11)         Mrl  ^  ^°~\   ■  -2a(t)Aj  n  do". 
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The  proportionality  factor  2o(t)  depends  upon  local  properties  of  the  surface, 
the  media  and  the  field,  in  accordance  with  the  localized  nature  of  the  radiation 
process.  We  will  consider  the  determination  of  a(t)  later. 
Now  we  integrate  (11)  from  t  =  0  and  obtain 

^  o 

The  ratio  do-  (0)/do-(t)  is  the  ratio  of  strip  widths  at  0  and  t,  or  more  precisely, 
the  limit  of  this  ratio  as  the  strip  width  tends  to  zero.  The  point  t  =  0  is 
the  point  of  tangency  Q,  so  A,(0)  =  A,(Q- ).  This  initial  value  of  A^  is  given  by 
(10),  We  now  insert  this  into  (12),  include  the  phase  given  by  (3)  and  denote 
dar(O)  by  do-(Q^),  Then  (12)  becomes 

ii/2        r   ^     t   ,   t 


(13)    u^(t)  =  D(Q-l)A^(Q^) 


n(Q^)  do- (Q^) 
^irn~   do-(t; 


exp  \ik 


K^S.^*  f   "^*|  ■   f"^*^'^*! 


Equation  (13)  gives  the  field  on  a  diffracted  ray  on  the  surface. 

Let  us  now  determine  the  diffracted  field  u,(P)  at  a  point  P  off  the  surface. 
First  we  denote  by  P,  the  point  on  the  surface  at  which  the  diffracted  ray 
through  P  leaves  the  surface.  Then  (13)  gives  u,(P, ).  Next  we  attempt  to  employ 
(6)  to  calculate  the  field  at  P  in  terms  of  that  at  P^.  This  cannot  be  done 
directly  since  the  surface  is  a  caustic  of  the  diffracted  rays  and  therefore  in  (6), 
do"  is  zero  on  the  surface.  Therefore^ instead, we  use  (6)  to  relate  vi^(P)  to 
u,(Q)  where  Q  is  any  point  near  P,  on  the  diffracted  ray  from  P,  to  P.  Ihen  we 

let  Q  approach  P, .  Since  dcry^do-  tends  to  zeio,  it  follows  from  (6)  that  A^s  ^d^*^^ 

1/2 
must  become  infinite  in  such  a  way  that  A  (do- /dor)  '  has  a  limit. 

If  3  denotes  the  distance  from  Q  to  P,  then  dcr^  =  dor  (Q)  tends  to  zero  like 

s.  Therefore  we  may  define  A^(P^)  and  da»(P,  )/da(P)  by 
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ilk) 


A,(P,) 


lira  s-"-/^  A,(Q) 


(15) 


dcr'{?j) /da- {?)  '     lim   s"'''  dcr  (Q)/dcr  (P). 


With  these  definitions  (6)  becomes 


(16)     u^(P)  =  A^(Pj^) 


n(Pj^)  dcr'(P^) 
lT(Tr  dcr(P) 


Tl/2 


exp  ik  ^i^d^^l^*  I  "t(s)J<^s 


Since  the  diffracting  surface  is  a  caustic  of  the  diffracted  rays,  dcr(P) 
is  zero  when  P  is  on  the  surface.  Consequently  (16)  does  not  apply  on  this 
surface.  In  fact,  the  present  theory  always  fails  in  this  way  at  a  caustic  and 
requires  a  modification  which  will  be  described  later,  Therefore,we  are  unable 
to  identify  Aj(Pt)  in  (16)  with  the  amplitude  of  the  diffracted  field  on  the 
surface  at  P, ,  given  in  (13).  Instead  we  assume  that  A^(P^)  is  proportional  to 
this  amplitude 

rn(Q,)  do-(Q  )  l^/^    r  >      -, 

(17)    W "( V"( VM«i>  STCT  3Fna ]     ==5.[-J°(^H- 

-  o 

In  (17),  which  deteiniines  A,(P^),  the  proportionality  factor  D(P^)  is  again 

called  a  diffraction  coefficient.     We  assume  that  it  is  the  same  function  of  the 

properties  of  the  field  and  object  at  P^  as  is  D(Q^)  of  the  properties  at  Q^. 

This  assumption  is  based  upon  the  reciprocity  principle  -  that  a  source  at  Q 

produces  the  same  field  at  P  as  a  source  at  P  produces  at  0.  Our  assumption 

assures  that  this  principle  is  satisfied  by  the  diffracted  field.  Now  we 

insert  (1?)  into  (l6)  to  obtain 

-Tl/2 
(18)   u^(P)  -  D(P-j^)D(Q^)u^(Q^) 


n(Q^)  d<y-(Q^)  dcr^P^^) 
liaT"  clo-(t'^)  dFTFT 


exp  I  (ilcn-<i)dt 

s 
f  ndtl  • 


+  ik 
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To  obtain  this  result,  we  have  necessarily  excluded  cafees  in  which  the 
field  or  some  component  of  it  has  to  vanish  on  the  boundary.  This  was  required 
because  our  deductions  were  baaed  on  the  amplitude  function  A ,(t)  of  the  field 
on  the  surface.  As  previously  noted,  the  surface  is  a  caustic  of  the  diffracted 
field.  Therefore  the  field  must  be  much  stronger  within  a  surface  layer  a  few 
wavelengths  thick  than  at  points  a  great  number  of  wavelengths  from  the  caustic. 
Thus  the  argument  used  to  derive  (18)  still  applies  if  A,(t)  is  any  measure 
of  the  amplitude  in  this  layer.  However,  the  argument  will  now  cover  the  cases 
>riaich  were  previously  excluded.  Within  the  layer,  the  field  will  be  characterized 
by  some  definite  variation  with  distance  from  the  surface,  At  any  point  other 
than  a  zero  of  this  profile,  the  field  amplitude  can  serve  as  the  measure  of 
the  field  in  this  layer. 

The  result  in  (18)  does  not  account  for  the  fact  that,  generally,  the 
field  within  the  caustic  layer  will  be  canposed  of  different  modes,  each  with 
its  own  profile.  Each  mode  will  then  have  its  respective  amplitude  A^  and 
thus  its  own  diffraction  coefficient  D  (Q^ )  and  proportionality  factor  o  (t). 
Therefore,  (18)  must  be  replaced  by 

-11/2         ^        t  s 


(19)        u^(P)  »  u^(Qi) 


Tx{Q^)       dcr(Q^)       dcr'(P^) 
'KCfT  '  do-  {?^}   *  d5=Tn" 


exp  Uc  j  ndt+ik  |    nds 


o  o 


r  1 

m  ^      L 


The  diffraction  coefficients  D  and  the  decay  exponents  a  must  now  be 

m  m 

determined.  These  coefficients  have  previously  been  studied  in  cases  of  diffrac- 
tion by  convex  objects  in  homogeneous  media'- •'''■  •'.  In  those  cases  it  was  found 
that  the  only  geometrical  property  of  the  surface  upon  which  D  and  a  depend 
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is  the  noimal  curvature  of  the  surface  in  a  plane  tangent  to  the  diffracted  ray. 

However,  in  those  cases  the  rays  not  lying  on  the  surface  are  straight  lines 

whereas  in  the  present  problems  they  are  curved.  The  physical  nature  of  the 

diffraction  process  suggests  that  in  the  present  case  D  and  a  will  depend  only 

upon  the  relative  curvature  of  the  incident  or  diffracted  ray  and  the  surface. 

3y  'relative  curvature'  we  mean  the  sum  of  the  curvatures  of  the  ray  and  the 

normal  section  of  the  surface  in  a  plane  containing  the  incident  or  diffracted 

ray.  Furtheimore  we  assvime  that  D  and  a  are  the  same  functions  of  this 

mm 

'relative  curvature'   as  they  are  of  the  surface  curvature  in  the  homogeneous  case, 

Of  course  the  value  of  the  index  of  refraction  n  which  appears  in  the  expressions 

for  D     and  a     must  be  that  at  the  point  at  which  D     or  a     is  to  be  evaluated, 
mm  mm 

Let  us  consider  a  scalar  field  satisfying  the  reduced  wave  equation  in 

an  inhomogeneous  medium  and  an  impedance  boundary  condition  on  the  body.     From 

the  known  values  of  D     and  a     for  such  a  field  in  a  homogeneous  medivira  these 

mm  *^ 

coefficients  can  be  evaluated  for  an  inhomogeneous  medium  by  means  of  the 
foregoing  assumption.    In  [2],[li]  the  coefficients  were  found  for  a  homogeneous 
medium  in  which  n  *  1.     To  obtain  the  corresponding  results  for  a  homogeneous 
medium  with  any  other  value  of  n  it  is  merely  necessary  to  replace  k  by  kn  in 
those  results.      Then  from  the  above  assumptions  we  obtain 

n3A  [kn(P^)]-^/^^[K(P^)-Kb-^(P^)]-^/^exp||g-  2i(3-\)^^^' 


(20)     D  (P.,) 


"  '      6VVA[v^(.v) .  ^.h.J]^'^ 


(21)      a(r^)  .  .-^"/^ 
m    1 


<^[6-Wp,f /3  |^K(P^)*b-l(P,)]2/3. 


-1/ 


In  these  equations  K(P-)  and  b~  (P,)  are  respectively  the  curvature  of  the  ray 
and  of  the  normal  section  of  the  surface  at  P,  .     ^i\)  is  the  i'iry  function  given 
by  (53),   Part  II  and  q^  is  the  m-th  zero  of   (^2),   Part  II. 
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U,  Plane  waves  in  a  plane  stratified  medium 

As  a  first  application  of  the  general  method  just  described,  we  consider 

an  unbounded,  plane-stratified  medium.  The  index  of  refraction  is  n  =  n(x)  and 

tends  to  the  constants  n  ,  n.  respectively  as  x  approaches  ioo.  We  suppose  that 

the  incident  field  at  x  =  -oo  is 

ikn  (x  cos  a  +  y  sin  a) 
(22)     u^  =  e   ° 


For  this  medium,  (l)  reduces  to 


/,-x     dV  _  n     d^x  ^  d  /  1  2n 
(23)     -i  -  0»    -2  =  dH^  ^  °  ^' 
dt         at 


A  first  quadrature  of  (22)  yields 


By  choosing  the  positive  sign  and  letting  x  tend  to  -oo,  we  deduce  that 

a  =  n  sin  a.  The  multivaluedness  of  the  slope  suggests  that  more  than  one  ray 

0 

can  pass  through  any  point.  The  number  of  rays  depends  on  the  nuiiiber  of  zeros 
X  of  the  equation 

(25)        n2(x^)  -  a^  -  0. 

We  call  points  at  which  (25)  is  satisfied  turning  points  since  at  them  the  rays 

become  vertical  and  are  turned  back  toward  -oo. 

We  will  consider  three  cases.   If  n  sin  a  is  less  than  the  minimum 

o 

of  n(x),  there  are  no  turning  points.  The  rays  then  form  a  family  of  congruent 
curves  proceeding  always  to  the  right.  If  n  sin  a  lies  between  the  miniimim  and 
maximum  of  n(x),  then  (25)  has  at  least  one  solution.  If  there  is  only  one 
solution  (turning  point)  x  =  x  ,  all  the  rays  become  vertical  along  the  line  x=x^ 
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and  are  refracted  back  toward  x  »  -oo.  The  region  x  >  x  is  then  a  refraction 
shadow  and  is  void  of  real  rays.  I'fe  will  also  consider  the  case  when  there  is  an 
additional  turning  point  at  x  =  x '  >  3C  • 

(a)  No  turning  points 

In  this  case,  we  inmediately  obtain  the  equations  of  all  the  rays  by- 
integrating  (2U),  choosing  the  positive  sign.  We  find 


(26) 


(  a  dx 

o 


If  we  use  the  value  of  a  in  terras  of  a,  this  can  be  put  into  the  foim 


(27)    y  =  X  tan  a  +  J 


n  sm  a 

_  (n^-n^  sin^a)^/^  "  ^ '^ 
-00  L    o 


dx  +  b. 


Thus,  for  large  negative  x,  the  rays  behave  like  straight  lines. 

To  find  the  phase  along  these  rays,  we  first  transform  (3)  to  the  parameter  t 
using  (2)  to  show  that  ds/dt  =  n.  Then  we  use  the  value  of  dx/dt  obtained  by 
integrating  the  second  of  (23),  and  we  find 


(28) 


0  '^^*   a(y-y^)  >  |  {nW)^^^dx. 


^o 


The  constant  ^    can  be  determined  by  comparing  (28)  with  u.  at  x  =  -oo.  Then 
(28)  can  be  written  in  the  form 


X 


(29)       0  ■  n  (x  cos  a  +  y  sin  a)  +  I  g(x)dx. 
Here,  g(x)  is  defined  by 


-00 
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2     2        2     l/2 
(30)  g(x)  -  (n  -n^sin  a)  '     -  n^cos  a. 

Now  we  find  the  amplitude  by  applying  (5).     By  noting  that  dcr  /dy=[l+(^)  J         > 
we  find  that 


(31) 


dcr^         dcr 


o_   /  dcr     _      n 
dy  /  dy      '  n(x  ) 


r     2  2       2" 

n  (x^)  -  n^sin  a 

2     2   .   2 

n  -n^sm  a 

0 


1/2 


Upon  using  (31)  in  (5)  and  letting  x^  ->  -oo,  we  obtain 

r  n2(x J-n^sin^a  1^^      A  (n  .os  a)^/^ 


(32) 


A  =  A(x^) 


o       o 

""5 — 5 7~ 

n  -n  sin  a 
o 


o^   o 


(n^-n^sln^c)^^  ' 


Comparison  of  (32)  with  u.  at  -co  shows  that  A^  »  1.  Finally,  by  combining  (28) 
and  (32)  we  obtain  for  the  field 

1/2  * 

(n  cos  o)  '  r  n  /^      /   ^    7 

(33)  u  =     -S— 5 TTTT  e^cp  iknjx  cos  a  +  y  sin  a)  +  Uc         g(x)dxl 

(n  -n^sin  a)  '^         ^  l^ 

This  result  agrees  exactly  with  (6)  of  Part  II,  which  was  obtained  by  asymptotically 
expanding  the  solution  of  the  corresponding  boundary  value  problem. 


(b)  One  turning  point 

Let  the  turning  point  be  at  x  =  x  .  The  incident  rays  therefore  become 
vertical  at  x  =  x  and  the  line  x  =  xl  is  then  a  caustic,  ^t   the  plane  x  =  x^, 

the  incident  ray  splits  into  a  real  ray  refracted  back  towards  x  =  -00  and  a 

2   2  1/2 

complex  ray  (see  Figure  2).  This  corresponds  to  (n  -  a  )  '  being  pure  imaginary 

for  X  >  X  so  that  there  are  no  real  solutions  of  the  ray  equation.  The  region 
a 

X  >  X  is  a  refraction  shadow  in  which  there  are  only  complex  or  imaginary  rays, 
a 

^   complex  or  imaginary  ray  is  a  complex-valued  solution  of  the  ray  equation.  We 
say  that  it  lies  in  a  given  region  if  the  real  part  of  x  lies  in  that  region. 
The  field  'along'  such  a  ray  is  found  by  utilizing  the  analytic  continuations  of 
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^x(r) 


(a) 


(b) 


Figure  1.  Typical  graphs  of  the  index  of  refraction  n(x)  or  n(rj  in  a  plane  or 
cylindrically  stratified  medium,  n  tends  to  n^  at  infinity.  If  the  index  is  of 
the  form  in  (b)  the  medium  will  contain  a  ductj  if  it  is  of  the  form  in  (a)  there 
will  be  no  duct. 


Im  X  ' 


Rex 


(b) 

Figure  2,  (a)  The  ray  pattern  in  a  plane  stratified  medium  due  to  a  plane  wave 
incident  from  the  left.  There  is  one  turning  point  at  x^.  The  incident  rays  have 
the  slope  tan  a  at  -oo.  (b)  The  ray  pattern  in  the  complex  x-plane.  The  real  r?y3 
become  vertical  at  x  =  x  and  then  turn  back  to  the  left.  The  region  x  >  x^^  is  a 
refraction  shadow  and  only  imaginary  or  complex  rays  occur  here.  For  Re  x  <  x^, 
thc?re  are  real  rays. 
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(3)-(6).  The  equation  of  a  complex  ray  is  the  analytic  continuation  of  the 
incident  ray. 

The  incident  rays  and  the  phase  and  amplitude  associated  with  them  are  again 
given  by  (26),  (28)  and  (32)  respectively.  In  each,  however,  we  replace  (x^jy^^) 
by  (x  ,y  )  and  all  three  formulas  hold  only  for  x  <  x„. 

The  equation  of  a  refracted  ray  after  turning  is  obtained  by  reflecting  the 
incident  ray  across  the  horizontal  line  y  =  y  through  the  turning  point  (3!^jy^)» 
The  ray  is  given  by  (26)  with  a  minus  sign  in  front  of  the  integral  and  i^tYg^) 
instead  of  (x  ,y  ).  This  result  can  be  put  into  a  foim  similar  to  (2?).  This 
form  again  shows  that  a  refracted  ray  behaves,  as  x  tends  to  -co,  like  the 
straight  line  with  slope  -tan  a. 

The  phase  along  a  refracted  ray  after  turning  is  obtained  by  using  (3)  with 
0  the  phase  of  the  incident  ray  at  x  =  x  .  We  thus  obtain 

x„ 

x  a 

(3U)    0  =  n^(-x  cos  a  +  y  sin  o)  -  1  g(x)dx  +  2(n^x^cos  a  +  I  g(x)dx)  . 

-00  -00 

We  recall  that  x  =  x  is  a  caustic  for  the  real  rays.  When  a  ray  (in  a  homo- 
geneous medium)  passes  through  a  caustic,  the  phase  of  the  field  is  retarded 
by  n/2  '^  ■' .   The  field  along  an  incident  ray  here  has  the  same  local  behavior 
at  X  «  X  .  With  this  phase  change  taken  into  account,  the  amplitude  along  a  refracted 

ray  is  given  by 

(n  cos  a)l/2e-"i/2 
(35)    A  -   ° 


(n^-n^sin^a)^/^ 


The  above  results, when  collected, yield  the  field  at  any  point. 
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(n^cos  a)^/2 
(n  -n  sin  a)  '^ 


ikn  (x  cos  a  +  y  sin  a)  +  ik  1  g(x)dx 

-00 


V  exp 


X  a 

ikn  (-X  cos  a+y  sin  a)-ik  I  s(x)dx+2ik  j  g(x)dx+21kn  x^coa  a  -  in  A 


-00 


-00 


»  x<: 


This  result  agrees  exactly  with  (11),  Part  II  obtained  by  asymptotically 
expanding  the  solution  of  the  corresponding  boundary  value  problem. 

In  the  refraction  shadow,  the  equation  of  a  ray,  the  phase  and  the 
amplitude  are  given  by  the  analytical  continuations  of  the  corresponding 
formulas  (26),  (28)  and  (32)  for  an  incident  ray.  The  field  at  any  point 
(x,y)  in  this  region  is  given  by 


(n^cos  a)^/2 
(n  sin  a-n  )  ' 


a 
ikn  y  sin  a-k  (x)+ikj   g(x)dx+ikn  x  cos  o-inA 


-co 


X  > 


This  agrees  with  (12),  Part  n,  the  asymptotic  expansion  of  the  exact  solution 
to  the  corresponding  boundary  value  problem, 
(c)  Two  turning  points 

Let  the  two  turning  points  be  at  x  and  x',  ^c  <  x'.  The  lines  x  =  x  , 
X  =  x'  are  caustics  for  the  real  rays  in  the  regions  x  <  x„,  x  >  x '  respectively. 

Q,  —   CI     —   cx 

We  will  label  these  half-spaces  I  and  III,  The  region  x  >  x^  is  a  refraction 
sliadow  and  in  the  region  between  the  two  caustics  (which  we  label  II)  there  are 
no  real  rays.  However,  we  will  nox^r  show  that  through  every  point  in  the  space 
there  pass  an  infinite  number  of  rays.  To  see  this  we  note,  that  at  the  plane 
X  =  X  the  incident  ray  splits  into  a  ray  infracted  back  into  I  and  an  imaginary 
or  complex  ray  continuing  into  II,  ;^t  the  plane  x  =  x',  this  imaginary  ray  splits 
into  an  imaginary  ray  refracted  back  and  a  real  ray  continuing  into  III,  This 
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splitting  occurs  each  time  a  complex  ray  reaches  either  of  the  planes  x  =  x 

a 

or  X  =  x^.  Each  of  the  additional  real  rays  which  arise  because  of  the  splitting 
of  a  cOTiplex  ray  is  coincident  with  one  of  the  initial  real  refracted  or  trans- 
mitted rays.  Consequently  there  will  be  infinitely  many  real  rays  passing  through 
each  point  in  I  and  III  and  infinitely  many  complex  rays  through  each  point  in  II 
(see  Figure  3*U). 

In  I,  the  equations  of  an  Incident  ray  and  its  continuation  beyond  the 
turning  point  are  the  same  as  those  equations  given  in  the  section  on  the  one- 
turning  point  case  and  in  the  region  x  <  x  .  The  equations  of  the  complex  rays 
in  II  are  the  appropriate  analytic  continuations.  The  rays  in  III  are  given  by 
(26)  with  (y-QsV^)   replaced  by  (xj,y^). 

Now  let  0  (P)  and  <ji  (.F)   stand  for  the  phase  at  a  point  P  along  a  ray.  The 
superscripts  t  and  r  mean  that  the  ray  at  that  point  is  being  transmitted  to  the 
right  or  to  the  left  respectively.  The  subscilpts  I,  n,  or  III  refer  to  the 
particular  region  in  which  P  lies.  The  subscript  n  indicates  the  number  of  times 
any  one  of  the  complex  rays  described  above (which  leads  analytically  to  a  ray 
being  considered) has  'reflected'  at  x  =  x^  and  x  =  x'.  Then  j0q  and  ^0^   are  given 
respectively  by  (29)  and  (3U).  j^q   is  given  by  the  analytic  continuation  of 


Ji^  ,  In  general,  we  have  for  m  =  1,  2,  ... 


X  ^a 


(37) 


]^2m  ~  "^o^"^  ^^"  °  "*"  ^a^°^  tt)  -  I   (n  -n^sin  a)  '  dx  +  I   g(x)dx 


X  -00 

a 


^c; 


d  I   (%sii 


+  2ini  I   (nSin^a  -  n^)^/^  dx. 


(38) 


X  ^     ^a 

^2m  "  ^o^y  ^^"  ^  *  ^a°°^  a)  +  i  J  (n^sin^a  -  x?)^''^   dx  +  j  g(x)dx 


Whi 


\  -«> 


K 


+   2mi  j  (n^sin^a  -  n^)"""/^  dx 


^a 
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(         o        0  0  1/2  I 

(39)       li^an-l  "     "o^^  ^^  ^  *  ^a'^*^^  a)  +  i   I     (n^sin  a  -  n  )  '   dx  +   I     g(x)<±3C 


-00 


.,       ..    r°-f  2   .   2  2sl/2    , 

Ksm-l)    I      (n  sin  a  -  n  )  '     dx 


^a 


(I4D)     ju02ni-l  "     ^o^^  ^^"  °  *  ^°°^  a)  +  J     (n^-n^sin^a)-'-4^  dx  +  J     g(x)dx 


^d  -°° 


a 
i(2m-l)   (    (n^sin^a  -  n^)^/^  dx. 


^a 


We  now  use  (5)  again  to  construct  the  amplitude,  making  the  same  assumption 
as  for  the  one-turning  point  case  and  an  additional  assumption.     Namely,  when 
a  complex  ray  splits,   the  amplitude  along  the  resulting  complex  ray  is  initially 
multiplied  by  i  e     '     .     Therefore,  if  we  define  ^   ,  A^  in  a  similar  way,  we  find 
that  jA*  and  jA^  are  given  respectively  by  (32)   and  (35)   and  that  jjA^  is  given 
by  the  analytic  continuation  of  (32).     In  general,  we  find  for  m  =  1,  2,   ... 

2m-l  ^^ 


"'^>     n^L-   ^■'"nA 


-2m       ,t     rani 
0  ® 


*To  justify  this,  we  consider  the  canonical  problem:   To  find  the  field  which 
is  outgoing  as  x  tends  to  -00  in  a  medium  with  an  index  of  Fip:ure  la  (with  only 
one  turning  point  at  x  =  x  ).  The  method  of  solution  parallels  that  for  the 
one-turning  point  case  and  jo-elds  the  amplitude  change  described  above. 
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21T1-1     ^ 

(^3)       n^2m-l  =     2  ijA^e 


/I  1  \  »t  o2-2m     .t     (m-l)rti 

^^^     III^2in-l=     2  A^ 


The  contributions  at  an  arbitrary  point  in  I,   II  or  III  form  geometric 
series  and  when  summed  yield  the  fields  in  the  respective  regions.     We  find 
that  in  region  I 


(n  cos  a) 


1/2 


7"2       2    .    2    xl/U 
(n  -  n^sin  a;   '^ 


exp 


ikn  (x  cos  a  +  y  sin  a)  +  ik   [     g(x)dx 


-00 


ikn  (-X  cos  a  +  y  sin  a)  -  ik  1     g(x)dx 

-oo 


x_ 


+  2ik 


I     g(x 

-CX) 


)dx  +  2ikn  x  cos  a  - 

o  a 


ITT 

T 


X  <  x„  . 

a 


Here,   ""^(x)   is  defined  by 

X 

ikS)       V|(x)     =     r     (n^sin^a  -  n^)-*-/^  dx. 


X 


In  the  interval  between  the  turning  points,   the  field  is 


(li7)       u 


1/2 
2(n^cos  a)  '     e 


kr^(x') 


(n^sui  a  -  n  ;  '    (1+ue  ; 


exp 


ikn  (y  sin  a  +  x  cos  a)+  ik 


a 
(    g(x)dxj 


-00 


exp[kiri(x)  *  ^]  +  2  exp^-k^Cx)  -^] 
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When  X  >  x' ,  we  find 
a 


ii(n  cos  a)  '     e 
o 


^^^^    "=  -7~1 — TIa"^ — 2k/i(x-)    ^^ 

o 


ikn  (y  sin  a  +  x  cos  a) 
0  a 


X  X 


+  ik  f  g(x)dx  +  ik  j  (n^-n^sin^a)-^/^  dx  . 


-00  X' 

a 

These  results  coincide  with  (11),  (13)  and  (ih)   of  Part  II.  The  latter  expres- 
sions are  obtained  by  expanding  the  exact  solutions  asymptotically. 

The  preceding  construction  has  been  done  in  a  heuristic  way.  By  looking 
at  the  complex-valued  solutions  of  (2U)  we  can  obtain  a  more  precise  description 
of  the  rays  and  associated  fields.  The  integrand  of  (26)  is  actually  a  double- 
valued  function  and  may  be  thought  of  as  being  defined  on  a  two-sheeted  Riemann 

surface  in  the  complex  x-plane  with  one  branch  point  occurring  at  x  =  x^.  This 

2  2 

corresponds  to  the  case  of  one  turning  point  whei-e  we  assvime  that  n  -a  has  a 

simple  zero  at  x  =  x  .  It  can  be  shown  that  the  integral  around  any  circle 

a 

in  the  x-plane  with  x  as  center  approaches  zero  as  the  radius  goes  to  zero.  Thus 
it  makes  sense  to  talk  about  integrating  up  to  x  and  then  from  x  (per>iaps  on 
the  other  sheet)  to  some  other  point.  Thus  y  is  a  double-valued  function  of  x 
ever^wl-ere  except  at  x  =  x  .  Since  only  real  values  of  x  are  of  interest,  the 
paths  of  integration  will  be  confined  to  the  real  axis  but  on  tv;o  sheets.  With 
this  interpretation,  (26)  with  (x  ,y  )  replaced  by  i^tJ^)   clearDy  gives  the 
two  branches  of  the  'real'  ray  corresponding  to  the  incident  and  refracted  ray 
for  X  <  X  .  That  is,  on  the  upper  sheet  we  take  (3U)  and  on  the  lower  slieet  we 
take  (3U)  with  a  negative  sign. 

The  Riemann  surface  associated  with  the  integrand  of  the  phase  on  a  ray 
is  the  same  as  that  just  described.  This  may  be  seen  by  substituting  in  (26) 
the  value  of  y  from  the  equation  of  the  ray,  Thei-efore^tbe  determination  of  the 
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radical  for  0  corresponds  to  the  determination  of  y. 

According  to  (32)  the  amplitude  is  a  four-valued  f\inctlon  of  x.     As  we  go 

2     2       2 

around  x     in  the  direction  in  which  arg  (n  -n  sin  a)  increases  we  can  follow  the 

determination  of  the  argvunent  of  A  modulo  n.     This  indetenninscy  is  resolved  by 
comparing  the  field  construction  with  the  asymptotic  expansion  of  the  exact 
solution. 

When  there  are  two  turning  points  at  x  ,  x',  the  Riemann  surfaces  of  the 
integrands  of  (26)  and   (26)  are  again  two-sheeted.     However,  because  of  the  two 
branch  points,  the  functions  y  and  0  are  infinite-valued.     At  any  point  P,   the 
phase  and  amplitude  will  be  given  by  (37)-(UO)   and  (Ul)-(l4li)   respectively.     A 
particular  value  of  m  mei«ly  corresponds  to  the  sheet  on  which  the  path  of 
integration  to  P  is  chosen. 


Figure  3 

The  ray  pattern  in  a  plane 

stratified  medium  due  to  a 

plane  wave  incident  from 

the  left.  The  incident 

rays  make  the  angle  a.   with 

the  X-axis  at  -co.     The 

rays  become  vertical  at 

x=x  and  then  turn  back 
a 

to  the  left.  Cnly 

imaginary  or  complex 

rays  occur  in  the  interval 

from  X  to  x",  and  at  xl 

they  give  rise  to  real 

transmitted  rays.  The 

values  X  and  xi  are  the 
a     a 

two  roots  of  the  equation 

2     2   2 

n  (x)-n  sin  aO  where  n 

is  the  value  of  n(x)  as 
X  tends  to  -co. 
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Rex 


Figure  h 
The  ray  pattern  in  the  complex  x- plane  in  a 
plane  stratified  medium  corresponding  to  the 
pattern  of  Figure  3.  The  lines  along  the  real 
axis  are  the  real  rays  and  the  curves  the 
complex  rays. 
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^.  The  field  of  a  line  source  in  a  plane  stratified  medium  with  a  plane  boundary 

In  this  section  we  construct  the  field  due  to  a  line  source  parallel  to 
the  z-axis  and  located  at  (x  ,0),  The  index  of  refraction  n  =  n(x)  tends  to  n^ 
as  X  approaches  infinity  and  satisfies  the  condition  n(x)  >  n(x^)  for  x  >  x^. 
The  medium  is  bounded  on  the  left  by  the  plane  x  =  x,  on  which  the  field 
satisfies  the  impedance  condition 

ih9)  ^  =  ikZu,      X  -  x^  <  x^. 

In  this  medium,  the  rays  are  determined  by  (2U).  The  constant  a  in 
(2U)  is  given  by  a  »  n(x  )  sin  a  where  tan  a  is  the  slope  of  a  ray  pt  the  source. 
Those  rays  for  which  0  <  a  <  n/2  proceed  to  the  right  without  turning.  They  are 
given  by  (26)  with  y  «  0.  Of  the  rays  for  which  ti/2  S  a  <  n  some  hit  the  boundary 
and  are  reflected,  while  others  are  turned  back  to  the  right  before  reaching  the 
boundary.  We  will  call  a  ray  from  the  source  a  direct  or  incident  ray  until  it 
is  refle  cted  or  tumedj  then  we  will  call  it  a  reflected  or  a  refracted  ray. 
For  one  particular  value  of  a,  say  a  =  a  ,  the  ray  grazes  the  bounding  plane. 
The  continuation  of  this  ray  is  the  boundary  of  a  refraction  shadow.  The  plane 
X  =  X,  is  the  other  boundary  of  the  shadow.  In  the  shadow  there  are  neither 
incident,  reflected  nor  refracted  rays  but  there  are  diffracted  rays  (see 
Figure   5 )•  All  the  rays  in  the  range  a  <  a  5  n  are  reflected  while  those 
in  the  range  n/2  <  a  <  a,  are  refracted. 

To  calculate  the  phase  along  the  incident  rays  we  apply  (3)  and  obtain 

the  result  (28)  with  y  and  0  both  zero.  To  determine  the  amplitude,  we  apply 

0      o 

(5)  and  obtain  an  expression  of  the  fonn  (32)  with  the  n\imerator  and  denominator 
respectively  multiplied  by  dy^  /da  and  By/Sa.  Calculating  these  derivatives  from 
the  equations  of  the  rays  and  denoting  by  A  all  the  constant  terms  that  arise. 
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we  find  that 


(51) 


.  .  A  fn-^-a^l-^A 


(   n  dx 

1    (^ 


a^ 


-1/2 


Here  x^  and  x^  are  respectively  the  larger  and  smaller  of  x  snd   x  .  The 
constant  A  depends  upon  the  strength  of  the  source.  If  the  line  source  is  of 
unit  strength  ,  then 


<«>   ^  '[A]  If^^o^  -  A 


t-iA 


Bj'-  combining  (5l)  and  (52)  with  the  phase  fonmula,  we  finally  obtain  for  the 
field  on  an  incident  ray 


(53)       u,  -  p(x,)-.^]-^^  [nW]-'^     8^  [   -2^ 

U  -I  '      (n  -a   ) 


1-1/2 


075 


exp 


ik   j     (n^-j 


a^)-'-^^dx  +  ikay  + 


in 
T 


Each  ray  for  which  o,    <  o  <  n  meets  the  boundary.     The  resulting 
i-eflected  ray  can  be  obtained  by  reflecting  the  incident  ray  across  the  horizontal  line 
y  =  constant  =  y.(:t )  through  the  point  of  reflection.     The  reflected  ray  is 
therefore  given  by 


By  a  line  source  of  unit  sti^ngth  at  (x  ,0),  we  mean  one  near  which  the  field 
behaves  like  the  Hankel  fvmction    iA  fU     [k:n(x  )r\  where  r  denotes  distance 
from  the  source. 
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($U)   y  -  2yj,(x^)  -  y^ 


fx 

u 

^ 

adx 

,  2  2n1/2  * 
(n  -a  )  ' 


^^ 


<  X  <  00. 


The  field  along  the  reflected  ray  can  now  be  obtained  by  using  (3)> 
($),  (8)  and  the  relation  u  (Q)  =  Ru.(Q).  The  value  of  do- /da  in  (5)  is 

r         1  O 

obtained  as  in  the  case  of  the  incident  rays,  3y  /3a  and  6y/6a  being  calculated 
from  {Sk)^     The  quantity  cos  5  can  be  determined  from  the  equation  of  the  incident 
ray  and  we  find  that 


(55)    cos  5  =  -  p(x^^)  -  a^  ^/^/n(Xj^) 


Combining  these  results  with  the  value  of  u.   at  the  point  of  reflection,  which 
can  be  obtained  froin  (53)  >  we  find 


(56)     u^  =  R   [n2(x^)-a2]-lA   [■„: 


2     2 


-lA 


Snk 


n  dx 


ITITpTf 


-1/2 


exp 


ilc 


(n^-a^)^/^^  +  ikay  +  ^ 


The  reflection  coefficient  R  is  given  by 


(57) 


[n2(x,)  -  af /^  >  Z 


TjITT 


[n  (x^,)   -  a  J 

The  expressions  in  (53)  and  (5?)  agree  exactly  with  (38)  and  (39)  of  Part  II, 
which  were  obtair^d  by  asymptotically  expanding  the  solution  of  the  correspond- 
ing boundary  value  problem. 

Those  rays  corresponding  to  the  parameter  range  "^  ^  o-   *^  %  turn  back 
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before  meeting  the  boundary.  The  equations  of  these  rays  and  the  field  along  them 

up  to  their  turning  points  are  determined  by  the  preceding  formulas  which  are 

now  valid  in  the  interval  x  <  x  <  x  .  We  note  that  (5l)  and  (^3)  are  indeterminate 

a  -   -  o 

at  X  •=  X  .  If  we  evaluate  the  limit  of  (5l)  as  x  tends  to  x  we  find  for  A  at 
a  a 


a  turning  point  the  re  stilt 


(58) 


iK(x  ) 

a 

8nk 


1/2 


p(x,)  .  a^] 


-lA 


Here  K{x  )  =  n'(x  )/n(x  )  is  the  curvature  of  the  ray  at  the  turning  point. 
Thus  the  field  at  a  turning  point  is  given  by  an  expression  of  the  form  (53) 
with  (58)  used  for  A  instead  of  (5l).  This  result  again  coincides  with  the 
asTTTiptotic  expansion  of  the  exact  solution  of  the  boundary  value  problem  given 
in  Part  II. 

Beyond  the  turning  point  of  a  refracted  ray,  the  equation  of  the  ray, 
is  given  hj   (5U)  with  x.  replaced  by  x  in  the  lower  limits  of  the  integrals  as 
we  see  from  the  symmetry  of  the  ray.  The  same  modification  yields  the  phase.  To 
find  the  amplitude,  we  use  (5).  As  before,  in  simplifying  this  expression,  we  have 
to  calculate  dy  /6a  and  ay/Sa  from  the  "equations  of  the  rays.  Differentiation  of  (5U) 
cannot  be  performed  directly  since  the  resulting  boundary  term  and  integral  are 
both  infinite.  One  way  in  which  this  difficulty  can  be  avoided  is  to  use  the 
smoothness  of  n(x)  to  integrate  (Sh)   by  parts.  The  resulting  expression  can 
then  be  differentiated  to  obtain  the  required  derivatives. 

Now  we  consider  the  shadow  i*egion.  The  diffracted  rays  which  cover  this 
region  are  produced  by  the  single  incident  ray  for  which  a  =  a^.  This  ray  is 
tangent  to  the  bounding  plane  at  the  point  (Xj^^Y^)  ■  Qj^.  The  resulting  surface 
ray  proceeds  along  the  boundary  x  =  x,  as  a  straight  line.  At  each  point  it 
sheds  a  diffracted  ray  so  all  the  diffracted  rays  of  (2U)  are  given  by  the  one 
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parameter  family  of  congnient  curves 


X 


(59)        y  =  I 


n(x,) 


dx 


[n^-nh^f/^    *    ^''' 


Here,   (Xj^jy^^)  =  P-,   is  the  point  of  shedding  of  the  rays  through  (x,y)  =  P. 

To  find  the  field  associated  with  these  rays,  we  now  utilize  (lU)-(21). 

The  surface  rays  are  parallel  lines  so  that  on  the  surface  d<r /dcr  =  1.  On  the 

o 

other  hand,  because  the  diffracted  rays  are  congnient,  we  find  on  evaluating 
the  limit  in  (l5),  that 

(60)    A^(P)  -  A^(P^)[n2-n2(x^,)]-^^  [n(P^)  K  (P^)]^/^ 

Since  the  diffracting  surface  is  a  plane  on  which  n(x)  is  constant,  the 
diffraction  coefficients  at  P,  and  Q,  are  equal  and  a  (f)  is  constant.  Further- 
more  all  other  quantities  in  (60)  may  be  evaluated  at  Q- ,  the  point  of  tangency 
of  the  limiting  ray,  instead  of  at  P, .  By  combining  (60)  with  (17),  we  find  for 
A^(P)  the  result 


(61)    A^(P)  -  A^(Q^)[n(Q^)  K  (Q^)]^/^  D2[n2-n2(x^)]-l>^  e^"^ 


b  ''b^ 


To  complete  the  calculation  of  u,  we  note  that  the  phase  change  along  the 
surface  ray  is  Just 

(62)    n(x^^)  (y^  -  Y^^). 

Furtheimore,  ^^(Qn)  can  be  determined  from  (58).  The  phase  along  the  diffracted 
rays  can  be  obtained  from  (5)  by  using  (62)  and  the  phase  of  the  incident  field 
at  Q- .  In  this  way  we  find 
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X     X. 


(63)    0(P)  =  n(x^)y  +  ]  f  ♦  f    [^^^-^^^^^j^^^  dx. 


If  for  any  point  in  the  refraction  shadow  y  <  0,  the  calculation  of  the  ray 
equation  through  P  and  the  field  at  P  can  be  carried  out  in  the  same  manner. 
Collecting  these  results,  we  finally  obtain  the  following  expression 
for  the  diffracted  field 


m 


n(x,) 


1/2 


exp 


f    X    ^^ 


ik  < 


The  diffraction  coefficient  D  and  the  decay  exponent  a  are  obtained  from  (20) 

m  m 


and  (21)  by  letting  b  ->  oo.  This  yields 


_  n3A[^(^)]-lA.  K  -'/\)e^\lt  -  2i(3-\) 


3/2 


(65)   D  = 


(66)   a  =  e 
m 


-/^  ^[6-Wx,)f  3,  2/3(^,. 


The  result  (6U)  with  these  values  for  D  and  a  agrees  exactly  with  the 
asymptotic  expansion  of  the  diffracted  field  obtained  from  the  exact  solution 
of  the  corresponding  boundary  value  problem  in  Part  II,  and  given  by  (56). 
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P(x,v) 


Refraction  Shodow 


P/XfV 


Figure  ^ 


Figure  $.     Diffracted  rays  in  a  re- 
fraction shadow  occurring  in  a  plane 
stratified  medium.     The  medium  is 
hounded  at  the  left  by  the  plane 
x:=x.    and  the  source  is  located  at 
(x  ,0).     Of  the  incident  rays  from 
the  source  some  are  reflected  from 
the  boundary  and  others  are  refracted 
back  to  the  right  before  reaching 
it.     One  ray  is   just  tangent  to  the 
boundary  and  the  region  above  it  is 
the  refraction  shadow.     This  region 
contains  no  incident  or  reflected 
ray  but  is  covered  by  diffracted  rays 
wMch  are  congruent  to  the  tangent 

^  ray.     They  originate 

from  the  tangent  ray 

at  its  point  of  tangency,   travel  along  the  boundary 

and  then  leave  it  tangentially. 


Figure  6.     Diffracted  rays  in  a  refraction  shadow 
occundng  in  a  cylindrically  stratified  medium. 
The  medium  is  bounded  internally  by  the 
cylinder  r=b  and  the  source  is  located  at 
(r  ,0).     Of  the  incident  rays   from  the 
source  some  are  reflected  from  the  boundary 
and  others  are  refracted  back  to 
the  right  before  reaching  it.     One 
ray  is   just  tangent  to  the 
boundary  and  the  region  above 
it  is  the  refraction  shadow. 
This  region  contains  no  incident 
or  reflected  ray  but  is  covered 
by  diffracted  rays  vrhich  are 
congruent  to  the 
tangent  ray.     They 
originate  from  the  tangent  ray  at  its  point 
of  tangency,   encircle  the  boumdaiy  any 
number  of  times  and  then  leave  it  tangentially. 


Figure  6 
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6.  The  field  due  to  a  point  source  in  a  plane  stratified  medium  with  a 
plane  boundary 

We  now  consider  the  field  of  a  point  source  situated  at  (x^,0,0).  The 
index  of  refraction  of  the  plane  stratified  medium  is  again  n  =  n(x)  and  the 
medium  is  bounded  on  the  left  by  the  plane  x  =  x,  on  which  the  field  satisfies 
(U9). 

The  rays  are  still  plane  curves  lying  in  the  planes  /y  +  z  «  0,  They 
are  determined  by  (23)  plus  the  equation 

2 

(67)  i|  =  0,    z  -  bt  +  2  . 
dt 

Now,  the  field  must  be  rotationally  symmetric  about  the  x-axis.  Therefore, 
for  each  value  of  /  =  —  -  const.,  the  rays  in  the  plane  ^^y  +  z  «  0  are 
characterized  by  the  second  equation  of  (23)  and 

(68)  r  =  a  t. 

Here,  a  *  a(l  ♦  J?   )  '     and  r    >=  y    +  z  .     Therefore,  we  can  make  use  of  all  the 
results  of  the  pirevious  section  if  we  merely  replace  a  by  a  and  y  by  r.     The 
only  modification  occurs  in  the  amplitude  calc\ilation.     The  amplitude  for  any 
direct  ray  is  given  by  (5)  multiplied  by  (r^/r)  '    .     do~^/do- is  still  interpreted 
as  the  limit  of  the  two-dimensional  area  ratio.     For  the  amplitude  we  thus 
obtain  (5l)  multiplied  by  r    '     and  with  a  replaced  byX     If  we  take  the  point 
source  to  be  of  unit  strength  ,  then  A     has  the  value 


By  a  point  source  of  unit  strength  at  (x  ,0,0),  we  mean  a  source  >riiich  produces 

ikn(x  )R 
locally  the  field    e  /UnR    where  R  is  the  distance  from  the  source. 
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A„  -  al/2(l«)-l   \fU^)  .  ~a2]-lA. 


Upon  taking  account  of  the  foregoing  remarks  we  finally  obtain  the  following 
expressions  for  the  incident  and  reflected  fields* 


(69) 


1/2 


^  =  ^  [nM"^^   P(x.)  .  .f-^ 


/ 


n  dx 


u   x< 


^;^:?7^ 


-1/2 


exp 


1/2 


^ 


ik 


/<n^- 


2\l/2^  ., 

a   }   '    dx  +  nJkar 


(70)    u^  .  i|Il  [„2..t^^  p(-.)-t'^ 


/      X  \ 

'  O         X 


r  i 


hi 


2, 
n  dx 


(717)^ 


-1/2 


exp  ik  <    I     +1     y  in  ~i 


a^)-'-/^dx  +  ikar 


To  find  the  field  along  the  diffracted  rays,  we  again  utilize  (Hi) -(21). 
This  time  the  surface  rays  are  radial  lines  in  the  plane  3C  =  x.  with  the  point 
(x.  ,0,0)  as  center.  Let  (jCj^^r^)  ■  P^  and  (x^^jR^)  -  Qj^  be  respectively  the  points 
of  shedding  and  of  grazing  incidence.  Then,  on  the  plane  dcr/do'*  '^b'^*  ^^ 
we  use  (5)  modified  as  above,  we  find  instead  of  (60)  that  the  amplitude  of  the 
diffracted  field  is  given  by 


(71)    A^(P)  .  A^(Q^)D2[n2.n2(x,)]"   [ 


-lAWxj,)  K  (x^,)r^, 


1/2  -^m^V^b) 


Now  along  the  surface  ray  the  phase  change  is 


(72)    nix^Xr^^B^), 


29 


Combining  (72)  with  the  phase  on  the  diffracted  ray  off  the  surface  and  with 
u.(Q,)  which  can  be  determined  from  (69),  we  obtain  for  the  diffracted  field 


<"'      "d  =  Z  ^  [n^.n2(:c,)]-'-A[„^(.„)-„^(x,,]-lA  „(^,  k  (^! 


exp 


X    ^o 


D  and  a  are  again  given  by  (65)  and  (66)  which  are  obtained  from  (20)  and  (21) 
m     m 

by  letting  b  ->  oo.  The  result  (73)  for  the  diffracted  field  agrees  exactly 
with  (62)  in  Part  11,  which  is  the  asymptotic  expansion  of  the  exact  solution 
of  the  corresponding  boundary  value  pixDblem. 

7.  The  field  of  a  line  source  in  a  cylindrically  stratified  medium  with  a 
circular  boundary 

Let  us  now  consider  the  field  due  to  a  line  source  parallel  to  the 
z-axis  and  located  at  (r  ,0).  The  index  of  refraction  is  n  »  n(r)  and  tends 
monotonically  to  the  constant  n,  as  r  ->  oo.  The  field  satisfies  the  impedance 
condition  (k9)   with  du/Sx  replaced  by  3u/9r  on  the  circular  boundary  r  =  b. 

In  this  medium  (1)  and  (2)  reduce  to 


(.     d^r    /d9\2  _  d  fl  2) 


(7U) 


and 


,,  X  ^  d^e  .  5  dO  dr    ^ 
(b)  r  -^  +  2  ^  ^  -  0 

at 
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A  first  quadrature  of  (7Ub)  can  be  perfonred  and  yields 

^"^^^        HIE    "    ^  • 
r 

This  can  be  inserted  into  (7Ua)  and  the  resvilting  equation  can  be  integrated  to 

yield 

Ar.  0        o2  1/2 

(77)      g  .  i  (n^  .  Ij  )   . 

r 

Thus  from  (76)  and  (77)  we  have 
in^\  '31'  .  *  r  /  2  2  2n1/2 

Again,  a  may  be  identified  with  the  ray  slope  tan  a  at  the  source  and  then 
a  «  ±  r  n(r  )sin  a,  (see  Figure  6), 

In  the  previous  problems,  the  zeros  of  n  -a  (turning  points)  played 
an  essential  role.  The  analog  of  such  a  zero  is  now  any  solution  of 

(79)  rV(r)  -  a^  -  0. 

These  zeros  include  the  points  of  closest  approach  of  a  r^  to  the  origin,  as 
well  as  inflection  points.  Consequently  for  e«ich  of  the  two  limiting  rays  which 
meet  the  circular  boundary  at  grazing  incidence,  the  points  of  tangency  are 
zeros  of  (79).  If  +  a.  corresponds  to  the  limiting  rays  (79)  yields 

(80)  a  "  bn(b)  -  ir  n(r  )  sin  o,  . 

From  (78)  we  find  for  the  equations of  the  incident  rays 
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(81) 


e 


adr 
I  -a  ; 


Here  r^  and  r^  are  respectively  the  larger  and  smaller  of  r  and  r  .  To  find 
the  phase  along  these  rays,  we  apply  (3)  in  conjunction  with  (76)  -  (78)  and 
ascertain  that 


(82) 


0  =  ae 


/  <H 


2  1/2 

)    dr. 


The  amplitude  is  now  found  by  using  (5).  First  we  note  that  along  an 
incident  ray,  we  have 


do-     dr 
d5  ■  ""H? 


-1/2 


t^  ^  <g'r  • 


By  making  use  of  this  expression,  we  find  that 
(83) 


d%  n[r^n^(x^)  -  a^]^/^(^y8a) 

^     ■     n(x  )[rV-a2JV2(9e/3a)   * 


The  derivatives  in  (83)  can  be  found  froro  (81).     If  we  insert  (83)  into  (5) 
and  lump  all  constant  terms  together,  we  obtain  for  the  amplitude 

-1-1/2 


(8U) 


1\  =  A. 


[.W]-A  f  f 


l''< 


m  dr 

/   2  2     2>3/2 
(r  n  -a  )  ' 


The  constant  A  in  (81j.)  is  proportional  to  the  field  strength  at  the  source. 
If  the  source  is  of  unit  strength,  then 

Upon  collecting  the  above  results,  we  obtain  for  the  incident  field 
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(86) 


u. 


[2  2/      X     2I-IA   [22     2I-IA 
Ir^n  (r^)-a  J     '^  [r  n  -a  J     '"* 


8nk 


2 

m  dr 

(rV-a2)3/2 


-1/2 


•exp 


Uc 


/[.3. 


2-|l/^ 
=j    dr  +  ikae  ■•■  ^ 


Those  rays  which  meet  the  boundary  non-tangentially  give  rise  to  a 

family  of  reflected  rays  with  an  associated  reflected  field.  The  details  of 

the  construction  of  this  field  %rill  not  be  given  since  they  ar«  similar  to  those 

employed  in  the  two  problems  of  Sections  5  and  6. 

Let  us  now  consider  a  ray  for  which  (79)  has  a  solution  r  =  r  .  Such  a 

a 

ray  has  r    as  its  distance  of  closest  approach  to  the  origin  and  after  reaching 
this  point  it  turns  away  from  the  origin.     The  equations  of  such  a  ray  and  the 
phase  and  amplitude  onituptor»r     are  given  by  (8l),    (82)  and  (8U)  respec- 
tively»     However,   (8U)  becomes  indeterminate  at  r  =  r  .     In  the  limit,  as  r 
approaches  r  ,   (8U)  yields 

£1 


(89) 


1  +  r 


n'(r„) 
a 

n(r  ) 

a 


1/2 


•    A^[l  -  r^  K  (r^)] 


1/2 


K(r  )  is  again  the  curvature  of  the  ray  at  r  =  r  .     On  the  outgoing  portion  of 
the  ray  for  which      r  >  r_        the  ray  is  given  by 


(90) 


9  = 


a  dr 

./   2  2  ^2a/2 
rvr  n  -a  ;  ' 


The  associated  phase  is 
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^'o 


(91)     0  -  ^  *  ]  j  +  [ 


*    ay 


9   «2  1/fe 
(n  -  -^  )    dr. 

r 


The  amplitude  is  given  by  an  expression  of  the  form  (8U)  with  the  integral 
replaced  by  the  svun  of  two  integrals.  The  limits  are  from  r  to  r  and  r  to  r 

respectively. 

We  now  turn  to  the  evaluation  of  the  diffracted  field  in  the  refraction 
shadow  formed  by  the  circular  boundary  and  the  continuation  of  the  limiting  ray. 
The  initial  direction  of  the  diffracted  rays  is  that  of  the  surface  ray  on  the 
boundary.  Ihis  ray  is  Just  the  boundary  itself.  All  the  diffracted  rays  which 
emanate  from  it  are  congruent.  Any  two  can  be  made  to  coincide  by  a  rotation 
of  one  of  them.  When  we  use  (80)  in  (81)  we  find  for  the  equations  of  the 
diffracted  rays 


•  r  n  -b  n  (bjl 


Here  (b,9.  )  =  P,  is  the  point  of  shedding  of  the  ray  through  (r,9)  =  P. 
Through  each  point  P  in  the  shadow  there  will  pass  a  double  infinity  of  diffracted 
rays.  The  surface  ray  which  arises  at  the  point  (b,Y,  )  =  Q,  will  encircle  the 
cylinder  infinitely  many  times  in  the  counterclockwise  sense.  After  each 
encirclement,  it  will  shed  a  diffracted  ray  at  P^  and  this  ray  will  pass  through  P. 
Another  surface  ray  will  arise  at  (b,-T  )•  Q,  going  clockwise.  It  will  also  shed 
diffracted  rays  at  (b,9jj)»  P^  which  will  pass  through  P. 

The  field  along  these  rays  is  obtained  by  using  (lii)-(19).  We  first 
determine  the  field  due  to  the  rays  which  travel  counterclockwise.  The  phase 
along  one  of  the  diffracted  rays  which  arises  after  the  surface  ray  has  made  q 
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counterclockwise  circuits  of  tiw  boundary  is  given  by 

r 
r 

(93)    0  =  bn(b)©  +  2nqbn(b)  +  •( 


b   b 


2     bV^b) 
n 5 

r 


1/2 


dr. 


Because  the  surface  rays  are  parallel  curves,  do-  /do"  »  1,  Upon  carrying  out 
the  limit  in  (15),  we  find  for  the  amplitude  on  a  diffracted  ray  the  result 

(9U)    A^(P)  -  A^(P^)[rV-bV(b)]~^^  [bn(b)  [K(b)+b-^j  ] 

To  determine  A,(P, )  we  note  that  the  index  of  refraction  and  the  radius 
a  J. 

of  curvature  are  constant  along  the  boundary  so  that  D  and  a  are  constant.  Upon 

mm 

evaluating  A,(P-)  and  inserting  it  into  (9U)  we  find  for  the  amplitude  on  a  ray 
which  leaves  the  surface  after  q  circuits  the  result 


-1}   1^/2 


1-1 A 


(95)    A^(P)  -  A^(Q^)D^rbn(b)  [K(b)  t-  b"^j      rrV-bV(b)]     expr-aj^(  r  ♦2TTqb)'l 


The  value  of  ^^(Q-,)  can  be  obtained  from  (85)  and  (89). 

When  we  collect  the  above  results  we  find  for  the  diffracted  field 
corresponding  to  each  mode  a  series  of  terms,  one  for  each  value  of  q.  This 
is  a  geometric  series  which  can  be  summed  explicitly.  Upon  summing  the  series 
we  obtain  for  the  diffracted  field  u,  produced  by  one  of  the  tangent  incident 
rays  the  result 


exp 


ik 


r      ^0 

b       b 


(87ik)- 


^2      bV(b)]^/^,     ,  .i 
n    --^7—        ^^*Y 


z 

m 


T)^     l-expbn  i  rbkn(  b  )+iba  1  lexprikbn(  b  )'^^^  ^%-\  A\ 
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In  (96)  the  diffraction  coefficients  D  and  the  decay  constants  a  are  given 
by  (20)  and  (21).  In  an  analogous  way  we  find  the  diffracted  field  uT  due  to 
the  other  tangent  incident  ray  which  gives  rise  to  the  surface  rays  encircling 
the  cylinder  clockwise.  We  obtain  an  expression  of  the  fonn  (96)  with  the  last 
exponential  in  the  sum  replaced  by 

(97)    expmcbn(b)(2n-e)  -  aj^b(2Tt  -  26  +  9^^  -  T^)|. 

The  total  diffracted  field  u.  is  given  by  Uj[  "  ^j  +  ^j« 

The  above  result  for  the  diffracted  field  agrees  exactly  with  (77)  of 
Part  II  obtained  by  asymptotic  e3q>ansion  of  the  exact  solution  of  the  corres- 
ponding boundary  value  problem. 

8,  Trapped  modes  in  a  plane  duct 

We  now  consider  a  medium  in  which  the  index  of  refraction  varies  as 

t 

in  Figure  lb.  If  a  source  placed  at  infinity  between  x  and  x  produces  rays 

with  a  given  value  of  a,  these  rays  remain  trapped  between  the  lines  x=x  and 

t  I 

x=x  .  The  points  x  ,  x  are  turning  points  for  the  rays  and  these  lines  are 

then  caustics  for  them.  Two  such  rays  will  pass  through  any  point  for  which 

X  <  X  <  X. .  This  corresponds  to  the  fact  that  when  the  parameter  a  <  n_  _, 
a       £t  max 

p  ?  1  /o  • 

the  quantity  (n  -a  )  '  is  pure  imaginary  for  x  <  x  or  x  >  x  .  Thus,  (2li) 

will  have  real  solutions  only  in  the  interval  x  <  x  <  x^  . 

Let  X  be  the  zero  of  n'(x)  which  lies  between  the  two  turning  points. 
The  two  rajrs  passing  through  (x  jO)  and  proceeding  towards  y."  00,  are 
characterized  by  the  equations 
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t 


I   *2[q/2]  j      .   (-1)^1/ 


adjc 


(98) 


/  X. 


y  =  S 


X  \ 


/  .  2[q/2]  (   .  (-1)^*^1 


l^a 


a  dx 


Here  [q/2]  is  ths  greatest  integer  less  than  or  equal  to  q/2,  and  q  is  an 
integer  related  to  the  number  of  times  the  ray  has  gone  through  the  caustic 
and  turned  back  after  leaving  (x  ,0)  and  before  passing  through  P  =  (x,y). 
The  field  at  P  is  the  sum  of  the  fields  on  the  two  rays  through  P 

11^,  (P)        iki2fp(P) 

(99)  VL  -  A^(P)e   ■"   +  A2(P)e   ^   . 

In  the  strip  x  <  x  <  x  ,V(j/  is  double  valued  and  0   has  infinitely  many  values. 
However,  only  two  of  the  branches  of  (jf   are  essentially  distinct,  the  others 
differing  from  these  by  additive  constants.  The  amplitude  A(P)  also  has 
the  same  behavior.  Now  the  field  given  in  (99)  must  be  single-valued  at  P 
and  we  assume  that  this  is  true  of  each  term  separately.  Then  the  difference 
between  corresponding  branches  of  the  functions  0   and  A,  60   and  5  log  A 
respectively,  must  satisfy 

(100)  ke(/  -  i  6  log  A  =  2ran  . 

Here  m  is  an  integer. 

^  covering  surface  can  be  associated  with  ^consisting  of  two  sheets  joined 
along  the  lines  defining  the  caustics,  since  along  these  lines  V<0    is  single- 
valued.  The  determinations  70,  and  V^^   ^^^  associated  with  points  on  the  upper 
and  lower  sheets  respectively.  On  the  same  surface  7  log  A  is  also  single  valued. 

If  t  is  the  unit  tangent  along  a  curve  on  the  covering  surface,  and  if 
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arclength  s  •=  0   corresponds  to  a  point  P  ,  then  at  any  point  "^  ^  ±3  given  by 

p)  r 


(101) 


0(P)  -  <jfi^^)  +   [         V0  '  tds  +  m(x)  70  •  ^  ds. 


The  last  integral  is  taken  around  any  closed  path  starting  at  P,  Therefore 
the  changes  in  (/  and  log  A  can  be  written 

(102)  60  -  (p  ^  •  t  ds 

(103)  6  log  A  =  0)  V  log  /  •  f  ds. 

The  change  in  log  A  corresponds  to  the  number  of  tijnes  a  ray  goes  through 
a  caustic  in  the  direction  of  increasing  phase.  In  making  a  complete  circuit 
along  the  cowering  surface,  a  ray  goes  through  a  caustic  twice  so  that 
5(log  a)  =  -n.  Inserting  this  and  (3D2)  into  (100),  we  obtain 

(lOU)    k  m  V0  .  f  ds  -  n  =  2nm. 


(p  ^ 


To  evaluate  the  integral,  we  take  the  closed  path  determined  by  the  ray 
which  goes  along  the  upper  sheet  from  P  on  x  =  x  to  P.  on  x  =  x  and  then 
along  th3  lower  sheet  to  P«  and  finally  back  along  the  caustic  to  P  ,  Now 
along  the  caustic 

^o 

(105)  j       V0  .  if  ds     -     ary(P^)   -  jCPj)]    • 

^2 
If  we  use  (98)  to  simplify  this  and  to  evaluate  the  remaining  integral,  we  find 


a  a 


x' 


I- 


2^2)1/2  ^. 


^a 
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Upon  inserting  (106)  into  (IQU)  we  obtain  the  equation 
t 

(107)    kj  (n^-a^)^/^  dx  =  (ra  +  ^)  n  . 

This  equation  determines  the  possible  values  of  a  for  each  given  value  of  m. 
According  to  (  99  )  and  (  IDl)  the  phase  of  each  terra  in  u  is  kay  +  ... 
Thus  ka  is  the  propagation  constant  in  the  y  direction. 

In  part  II  the  problem  of  trapped  modes  is  formulated  as  an  eigenvalue 
problem,  the  eigenvalue  being  the  proportion  constant  in  the  y  direction  divided 
by  k.  The  asymptotic  solution  of  this  eigenvalue  is  found  to  be  given  by  (19) • 
This  equation  coincides  exactly  with  the  above  equation  for  a.  This  agreement 
constitutes  a  partial  justification  of  the  geometrical  method  introduced  here» 
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Part  II.  Asymptotic  Expansion  of  Solutions  of  Boundary  Value   Problems 
1,  Introduction 

In  Part  I,  a  geometrical  method  was  given  for  calculating  fields.  The 
method  was  also  used  to  determine  the  fields  explicitly  in  various  diffraction 
problems.  In  this  section,  we  consider  these  same  diffraction  problems  formulating 
them  as  boundary  value  problems  for  the  equation 

(1)        [/\+  kV(x)ju  =0. 

We  then  solve  these  problems  explicitly  and  expand  the  solutions  asymptotically 
with  respect  to  k  as  k  ->oo.  The  leading  tenn  of  the  asymptotic  expansion  which 
we  call  the  asymptotic  form  of  the  solution  is  all  we  consider,  in  general.  We 
find  in  every  case  that  the  asymptotic  form  of  the  solution  of  the  boundary  value 
problem  agrees  exactly  with  the  field  determined  by  the  geometric  theory.  This 
agreement  constitutes  a  partial  demonstration  of  the  validity  of  the  geometric 
method  of  Part  I, 

Ir\   calculating  the  exact  solution  asymptotically,  the  \<IKB   method  and 
tlie  method  of  stationary  phase  are  applied,  Wliere  the  solution  is  given  by 
an  integral,  the  method  of  stationary  phase  yields  the  asymptotic  form  of  the 
field  whenever  there  exists  a  stationaiy  point.  This  method  gives  the  field  in 
the  lit  region.  For  those  field  points  which  lie  in  the  refraction  shadow, 
there  are  no  stationary  points.  We  then  transform  the  integral  into  a  Fourier- 
Bessel  series  by  making  a  Watson  transformation.  The  subsequent  analysis  then 
yields  the  field  as  a  sum  of  modes. 

In  the  analysis  of  diffraction  in  inhomogeneous  media,  similar  techniques 
have  been  employed  by  C,  L,  Pekeris'-  -^  for  an  acoustical  point  source,  by  F,  G, 
Friedlander*-  •'  for  a  pulse  and  by  N,  A,  Haskell'-  ■».  D,C,  Pridmore  Brown  and  A, 
Ingard'--'  have  verified  that  the  asymptotic  form  agrees  with  the  measured 
diffraction  field  in  the  refraction  shadow  for  sound  in  a  stratified  medium 
with  a  constant  temperature  gradient. 
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2,  Plane  waves  and  trapped  modes  in  a  plane  stratified  meditan 

Let  us  first  consider  a  plane  stratified  medium  for  which  n  =  n(x).  We  will 
assime  that  n(x)  tends  to  the  constant  n  as  x  tends  to  -oo  and  to  n,  as  x  tends 
to  +00.  Now  suppose  a  wave  exp[ikn  (x  cosa  +  y  sinaj]  is  incident  from  -oo.  Then 
the  field  u  will  be  proportional  to  exp[ikn  y  sina]  and  independent  of  z.  There- 
fore we  write  u  in  the  form 


(2) 


u(x,y)  ■  e 


ikn  y  sina 

f  (x) 


We  suppose  that  at  -oo,  f(x)  behaves  like 


(3) 


f(x)~e       ° 


ikn     xcosa 


+  R  e 


-ikn    xcosa 
o 


for  some  value  of  the  constant  R.     At  x  =  +oo  we  require  f (x)  to  satisfy  the  radiation 
condition J 


(U) 


lii^lf     -  iko,   cosa  f| 


X-5> 


Upon  inserting  (2)  into  (1)  we  find  that  f  satisfies 


(5) 


f   +  k^(n^  -  n     3in^a)f  »  0  . 
XX     ^      o      ' 


We  now  seek  that  solution  of   (^)  which  satisfies   (3)  and   (U).     We  will 
obtain  the  asymptotic  form  for  large  k  of  that  solution.     Then  we  will  insert 
it  into   (2)  to  obtain  the  asymptotic  form  of  u(x,y). 

A  point  at  '.^ich  the  coefficient  of  f  in  (5)  vanishes  is  called  a  turniniL^ 

2         2         2 
point.     If  there  is  no     turning  point  and  n    >  n     sin  a,  then  the  asymptotic 

form  of  f  can  be  found  by  the  WP33  method.     Upon  inserting  it  into   (2)  we  or.t?in 


(6) 


u(x,y) 


,n     cosa V 1/2 
(  o  ) 

(n  -n     Sin  a) 


exp 


ikn  (x  cosa  +  y  sina^ 


ina)  ♦  ik  I 

'"CO 


g(x)dx 
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Here  g(x)  is  defined  by 


(7) 


f    \        /  2    2.2  xl/2 
;(x)  =  (n  -  n  sin  a)  '  -  n 


cosa 


2         2         2 
If  there  is  just  one  turning  point  at  x  =  x     and  if  n    >  n     sin  a  for 

X  <  x_  then  for  x  /  x    the  asymptoticr  form  of  any  solution  of    (5)  is 


(8)     f(x)' 


(n^-n^  sin^af-^^  fk  cos(ki^(x)-f)   +  B  cos(k;i(x)+J)l 
(n^  sin^a-n^r^   |  e'^'^^-)   .  B  e^^^^^^ 


X  <   X 


X  >   X  • 

a 


Here  A  and  B  are  constants  while  ^  and '^ are  defined  by 
(9)  KM  ' 


J 


(n^-n^  sin^a)^/^  dx 


(10) 


1 


(x) 


/  2  .  2   2x1/2  , 
(n  sm  a-n  ;  '  dx  . 


Upon  imposing  the  radiation  condition  (U)  we  see  that  B  "  0.  Equation  (3) 
determines  the  value  of  A»  When  these  constants  are  inserted  into  (8)  and 
the  result  is  substituted  into  (2),  the  following  expansion  of  u  for  x  /  x  results 


(n     cosa) 
(n  -n     sm 


ikn  (xcoso  +  ysina)  +  ik  g(x)dx 

-00 


+  exp  ikn  (-xcosa+ysina)' 


X 

-ikl       g(x)dx+2ik| 

'-00  ^ 


)dx+2ik|       g(x)dx+2ikn  x  cosa-  in/U 

00 


f,      X  <  x^    . 
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(n  cosa)  '  r~ 

(12)  u-— -I ^  ,  A      expli 


)  u- — X—  exp  ikn    yslna  -  ki?{x)+2.k\       g(x)dx+ikn  x  cosa-inA  >      x    >  x  • 

^      2  .  2   2.lA        o         ;:'     J   °^       0  a  a 

n  sm  o-n  )  '^     L  ioo  J 


Now  let  us  consider  the  case  in  which  there  are  exactly  two  turning  points, 

one  at  x  and  the  other  at  x*  >  x  .  Then  the  second  asymptotic  form  of  f  in 
a  a    a 

(8)  applies  only  for  x  <  x'.  Two  other  forms,  corresponding  to  those  in  (8), 
apoly  in  the  regions  x'  <  x  and  x  <  x  <  x'  respectively.  By  matching  the 
constants  of  the  two  forms  which  both  apply  in  the  interval  between  the  turning 
points,  and  applying  (3)  and  (U),  we  completely  determine  the  asymptotic  form 
of  f.  When  this  is  inserted  into  (2),  the  asymptotic  form  of  u  is  obtained. 

In  the  region  x  <  x  this  form  can  be  obtained  from  (11)  by  multiolyin^  the 

""  -inAf   -2k^(3c«)1r  -2k^(x')']-l 

second  exponential  by  the  factor  e"  '  [U-e     °'  j[U+e     °  J  .  In  the 

interval  between  the  turning  points  we  have 


(13) 


u  -^ 


2(n  cosa) 
0 


1/2  ^*?(^a^ 


(n  sin  a  -  n  )  '   Ql+Ue   "-^  a  ) 


exp 


ikn  (ysina  +  x  cosa)+ik 


c 

-00 


g(x)dx 


•  <  exp 


rk.2(x)+inAl   +  2  exp   pk>^(x)-inAl  I 


For  X  >  X •    we  have 


(lU) 


-CO 


+ikf       (n^-n' 


2    .  2  a/2 , 

sm  a;       dx 


■] 
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These  results  have  been  ccwipared  with  those  of  the  geometrical  theory  in  Part  I, 
Section  U. 

Now  let  us  consider  the  trapped  modes  in  a  plane  stratified  medium.  These 
modes  are  solutions  of  the  form 

(15)  u(x,y)  -  e^^^y  f (x)  . 

In  (15)  a  is  a  constant  which  is  to  be  determined*  The  product  ka  is  called 
the  propagation  constant  of  the  mode  in  the  y  direction.  The  function  f (x)  is 
seen  to  satisfy 

(16)  f_  +  k^(n^-a^)  f  -  0. 

We  also  require  f (x)  to  satisfy  the  radiation  conditions 

(17)  3,11^00  l^x  -  ^^   (4-^^)^'^^  ^1  "  0 


lijn  _If^  +  ik  (n?-a^)^/2  f  |  -  0. 


^   '  X  -*>»00   X 

It  is  clear  that  this  problon  will  have  non-zero  solutions  only  for  certain 

values  of  a,  which  are  called  eigenvalues. 

2 

We  assume  that  n  (x)  has  the  behavior  shown  in  Figure  3b,  with  n  <  n, . 

It  is  readily  seen  from  the  asymptotic  forms  of  f (x)  constructed  above  that 

2    2     2    2 
asymptotically  there  is  no  solution  of  (16)  -  (18)  if  a  <  n^  or  a  >  n^g^^» 

Therefore  we  assume  that  n^  <  a  <  n"^  .  Then  there  ai-e  exactly  two  turning  points 

in  (16).  We  now  construct  the  asymptotic  form  of  f (x)  in  the  same  way  as  before. 

The  two  asymptotic  fomis  which  both  apply  between  the  two  turning  points  do  not 

2 
match  unless  a  satisfies  the  equation 

x' 

(19)  k(  ^  inW)^^^  6x'L  +^n   . 


-  lil;  - 


2  2 
Here  x  and  x'  denote  the  turning  points,  which  are  the  two  roots  of  n  -a  "0, 
a     a  ' 

and  m  denotes  any  non-negative  integer.  Equation  (19)  is  the  eigenvalue 

equation  and  its  roots  a  are  the  eigenvalues. 

The  eigenvalue  equation  (19)  can  be  solved  asymptotically  for  fixed  m  and 

k  large.  For  as  k  becomes  infinite  the  integral  must  tend  to  zero.  Since  the 

integrand  is  positive,  except  at  the  endpoints,  the  length  of  the  interval  of 

integration  must  tend  to  zero.  Therefore  x  and  x'  must  both  tend  to  x   « 

a  a  max' 

the  value  of  x  at  which  n  assiimes  its  maximum.     Upon  making  use  of  this  fact 
we  obtain  the  restilts 

(20)  a    -^  n        -  f'Jjik  y^  !!1-Li^ 

m        max      I       n  /  k 

V      max   / 


The  value  of  x     is  found  to  be 

(21)  X  ^x        -   (S!L^)^^Vn        n..   )"  ^^ 

a        max       ^       T^^       ^  "max    max^ 

The  value  of  x^  is  obtained  from  (21)  by  changing  the  sign  of  the  second  tem. 
The  asymptotic  foms  of  f  (x)  and  of  u(x,y)  are  now  easily  obtained  but  will 
not  be  given  here.  The  above  results  for  a^,  x^  and  x-  have  been  compared 
with  those  of  the  geometrical  theory  in  Part  I,  Section  8. 


^'     The  field  of  a  line  source  in  a  plane  stratified  medium  with  a  plane  boundary 

Let  us  again  consider  a  plane  stratified  medium  bounded  by  a  plane  at  x-x.  . 
Suppose  a  line  source  parallel  to  the  z-axis  is  located  at  (x  ,0).  Then  the 
field  u  will  be  independent  of  z.  As  before,  the  index  n(x)  tends  to  n^  as 
X  tends  to  +00.  The  field  again  satisfies  (1)  with  the  inhomogeneous  tern 
-  5(x  -  x^)  6(y)  on  the  right  side  for  x  >  x^^.  On  the  boundary  we  require 
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that  u  satisfy  the  impedance  condition 


(22)  ^'^'^Zu  . 


Here  the  constant  Z  is  the  impedance  of  the  boundary  plane.  If  p  denotes 
distance  from  the  source  in  the  x,y  plane  then  as  p  ->oo,  u  must  satisfy  the 
radiation  condition 

(23)  pll^^  P^/^  |u  -  ikUpl  -  0  . 

To  solve  (1)  (with  the  source  term),  (22)  and  (23)  for  u  we  introduce 
f(x,a),  the  rburier  transform  of  u(x,y)  with  respect  to  y,  defined  by 


CD 

(2k) 

'-00 


CD 

f(x,a)  -  j  u(x,y)  e"^^  dy 


From  the  formulation  of  the  problem  it  follows  that  u(x,y)  *  u(x,  -y).  Then 
(2U)  shows  that  f(x,a)  -  f(x,  -&)•  Upon  taking  the  Fourier  transform  of  (1) 
(With  the  source  term)  we  find  that  f (x,a)  satisfies  (16)  with  the  inhomogeneous 
term  -6(x-x  )  on  the  right  side.  Let  f^  and  f-  be  two  linearly  independent 
solutions  of  the  homogeneous  equation  (16)  which  are  even  functions  of  a*  We 
require  f,  to  satisfy  the  radiation  condition  (17)  and  f^  to  satisfy  the  boundary 
condition  (22).  Then  we  attempt  to  construct  f  in  the  form 


(25)  f(x,a)  -  i 


/c^  f^  (x,a)        ^o  -  ^  '^^^ 


^Cgfg  (x,a)        \<'^<\    , 


This  form  for  f  satisfies  (16)  for  x  /  x  ,  and  satisfies  (17)  and  (22). 

By  taking  account  of  the  term  -6(»-x  )  in  the  differential  equation  for 
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f,  and  requiring  f  to  be  continuous  at  x=3Co  we  can  determine  the  constants  c^ 
and  Cp  in  (25).  Then  (25)  becomes 

(26)  f  (x,a)  -  f^{y^,&)   fg^V^^  ^'^^^^       * 

Here  x  and  x  indicate  respectively  the  larger  and  smaller  of  the  values  x 
and  X  .  The  function  W(a)  is  the  Wronskian  of  :^(x,a)  and  f2(x,a)  which  is 
independent  of  x  but  which  depends  upon  a.  According  to  (26)  and  the  evenness 
of  f-  and  f,,  f  is  even  in  a  as  it  must  be.  Now  we  take  the  inverse  transform 
of  f (x,a)  and  obtain 

00 

f(x,a)  e^^^y  da  . 


(27)  u(x,y)  -  ^ 


00 

k 


/-oo 

It  may  be  verified  by  direct  computation  that  (2?)  is  the  solution  of  the 
problem.  For  our  subsequent  calculations  it  is  convenient  to  write  u  ■  I  +1 
where 

(28)  ^  •  ^  f  ^(^'^^^  ^^^^^  ^• 

'O 

We  win  now  evaluate  u(x,y)  asymptotically  for  large  k.  Two  different 
methods  of  evaluation  will  be  necessary,  one  for  points  x,y  in  the  illuminated 
region  and  the  other  for  points  in  the  shadov;  region.  For  points  in  the  illuminated 
region  the  integrals  in  (28)  will  be  evaluated  by  the  method  of  stationary  phase. 
The  stationary  points  will  be  found  to  correspond  to  the  rays  through  the  point. 
One  will  correspond  to  the  direct  ray  from  the  source  to  x,y  and  another  to  the 
reflected  ray.  For  points  in  the  shadow  region  there  will  be  no  stationary 
points.  Therefore  the  integrals  will  be  evaluated  differently,  namely,  in 
terms  of  residues  at  the  poles  of  the  integrands  in  the  complex  a-plane. 

To  perform  the  stationary  phase  calculation,  we  must  now  evaluate  f(x,a) 
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asymptotically  for  large  k.     Vfe  do  this  by  first  evaluating  f,   and  f^  asymptoti- 
cally as  we  did  in  the  previous  section.     For  this  purpose  it  is  convenient  to 
subdivide  the  a-axis  into  the  three  intervals  0  to  n(x,  ),  n(x,  )  to  n-   and  n-   to   oo. 
In  the  first  and  third  intervals  there   are  no  turning  points  for  any  value  of  x, 
while  in  the  second  interval  there  is  exactly  one  turning  point  x     for  each  x. 
Therefore,   for  a  in  the  first  interval 

X 

(29)  f3^(x,a)  ^  cCn^-a^r^^  expfik  (    {n^^^)^^^dx 

X 

(30)  f2(x,a)  -^  c'(n2-a2)"^^  sin[k  ^      {n^-a^)'^^^dx  +  y  . 

Here  c  and  c*   are  arbitrary  constants  and  T  is  a  solution  of  the  equation 

.1-1/2 


(31)  cot  Y    =     iZ   rn^(Xj^)-a^ 


For  a  in  the  interval  n(x,  )  <  a  <  n^,  there  are  two  cases.  If  the  turning 
point  x^  <  x^,  f^  is  given  bv  (29).  When  x^  <  x  <  x^,  (30)  and  (31)  still  apply 
for  f^  provided  n  -a  is  replaced  by  a  -n  ,  sin  replaced  by  sinh  in  (30)  and  cot 
replaced  by  coth  in  (31).  For  x_  <  x  <  x  ,  we  use  this  and  (8)  with  n  sin  a 
replaced  by  a  to  show  that 

/k/  (a2-n2)^/2dx  +  Y       ^ 
(32)   f^  ~  c'(n2-a2)-lA   e  "^  sin  ^k  J  in^^h^^^6^   .  j] 


X  ""a 

^a 


-k(  (a2-n2)l/2ax  -  Y       ^ 

w  e  cos 

'   -a 


(k  f  (n2-a2)l/2ax  *  j] 


On  the  other  hand,  if  x_  >  x^  >  x,  f,  is  again  given  by  (30)  and  (3l)  with  tte 

cl     O         c. 
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above  changes,  f,  is  given  by  (29)  for  x  >  x  provided  x.  is  replaced  by  x 
For  X  <  X  <  XL  J  we  again  use  (8)  and  shDw  that 


a 


(33)     fg  -  c(a2-n2)-l^ 


f     -k/  '(a2-n2)l/2dx  .  f   k/  '(a2.n2)l/2dx  -  f    ^ 
^  e  +  e  f    . 


Let  the  integrals  corresponding  to  each  of  the  a-intervals  be  denoted  by 
I*,  I*,  I*  respectively.  The  integrand  of  I*  may  be  expanded  asymptotically  by  using 
(29)  ani  (30)  for  f^  and  fg  in  (26),  This  yields 


(31''  ^-w;] 


X  , 
X   o 


J'^^Hexplikj   (n2-a2)^/^dx+ilcay)  *Rexp(ik  S  (  ^  /  (  (n^-a2)^/^dx+ikay^c 


o.?F 


[n^CxJ-a^l^^     I7Z^ 


Here,  the  quantity  R  is  defined  by 


(35) 


We  nov;  apply  the  method  of  stationary  phase  to  evaluate  each  of  the  two  integrals 
of  which  I,  is  the  sum.  The  stationary  points  a  are  given  respectively  by 


■Is 


(36)     "   •     ^  ^ 


11^^ 


=  0, 


(37) 


'   X   ^o 


a  dx 


,  2  „2>l/2 
(n  -a  )  ' 


0. 


If  the  point  (x,y)  is  in  the  lit  region  and  y  S  0,  these  two  conditions  can  be 
satisfied  for  certain  a.  This  is  apparent  because  (36)  and  (3?)  are  the  equations 
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of  the  incident  and  reflected  rays  discussed  in  Part  I,  Section  5.  If 

these  values  of  a  lie  in  the  first  interval  then,  as  we  will  see,  there  are  no 
stationary  points  in  I-,  I,  or  I~«  Therefore  I  is  given  by  the  sum  of  the  two 
contributions  from  the  two  stationary  points  of  I, ,  Denoting  these  contributions 


by  u.  and  u  we  have  u  «-•  u .  +  u  where 
1     r  1    r 

x>  '] 

ik  (     (n^-a^)^/^  dx  +  ikay  +  ^ 

^1       ^ 


exp 
(38)    u,  ^  


(l(x^-a2)^A  („2^2,lA/3„,  I  r^^ 


\         X^ 


(39)    Uy 


"XX 

R  exp  ikjj   +  r  I  (n^-a^)-'-/^dx  +  ikay  +  ^ 


The  asymptotic  form  of  1^,   where  the  turning  point  x  <  x  and  x  <  x,  is 

c  cL      O         cL 

given  by  (3U).     However   R   must  be  replaced  by  e"     ^     and  x.    by  x     in  the  lovrer 
limit  of  the  integrals  in  the  integrand  of  I^ .     With  these  changes,  the  stationary 

points  are  again  given  by  (36)  and  (37).     If  the  roots  of  the  equations  lie  in  the 

+  +  + 

second  interval  then  I^  has  two  stationary  points  while  L.  and  I-  have  none.  In 

this  case  the  asymptotic  form  of  u  is  again  given  by  (38)  and  (39). 

When  x  >  x  the  asymptotic  fonn  of  I^  is  obtained  by  using  (29)  and  (33)  for  f, 
and  the  modified  form  of  (30)  for  fp.  However,  on  applying  the  method  of  stationary 
phase  to  the  resulting  integral,  we  find  that  there  are  no  stationary  points.  Either 
(36)  and  (37)  correspond  to  parameter  values  a  for  which  there  are  no  rays  or  the 
integrand  contains  only  real  exponentials.  Finally,  the  asymptotic  form  of  I_  for 
any  (x,yj  also  involves  real  exponentials  so  that  we  obtain  no  field  contributions 


from  this  integral. 

When  y  <  0  the  saroe  results   (38)   and  (39)  are  obtained  from  stationary  points 
of  I",  while  I     has  none.     The  result  u    ~'  u.+u     has  been  compared  with  the  geometrical 
result  in  Part  I,  Section  5. 

Let  us  now  consider  a  point  (x,y) in  the  shadow  region.     For  such  a  point  none  of 
the  integrals  considered  above  has  a  stationary  point  so   the  method  of  stationary 

phase  does  not  apply.     In  order  to  obtain  a  representation  of  u  which  will  turn  out  to 

+ 
be  useful  in  the  shadow  region,  we  will  rewrite  I".     From  Fif^re    7    we  see  that 


(UO)  f  fe^^y  da  =         f    fe^^y  da  i  2T,i  Z.  K, 

OP  r+  0Q+  +  ^+ 


^  similar  formula  holds  for  the  integral  of  fe         "•  .     Here  B_^  denotes  the  possible 

+ 
branch  cuts  of  f(x,a)  in  the  half-plane  +Im  a  >  0  and  R~  denotes  the  residue  at  the 

mth  pole  of  the  integrand  in  the  same  half -plane.    In  the  ^pendix  it  is  shown  that 
for  ty  >  0  the  integral  of  fe"     ^^  over   P  tends  to  zero  as  the  radius  of  the  semi- 
circle becomes  infinite  while  for  ly  <  0  that  over   ^  "tends  to  zero.     In  the  Appendix 
it  is  shown  that  the  integral  of  fe"     ^^  over  0Q_  is  the  negative  of  that  of  fe       ^ 
over  OQ     asymptotically  for  large  k.      It  is  also  stown  that   the  branch  cut  integrals 
are  asymptotically  zero. 

Upon  making  use  of  the  above  results  in  (2?)  and  (28),  and  then  using  (26)  we 
find  that  u  is  asymptotically  given  by  the  residue  series 


(Ul)  u(x,y)   '~'  -ik  ^ 


m' 


m  ew(a^) 

ea 

In  (I4I)  the  a  are  the  roots  of  the  equation 

(h2)     W(a)  =  0. 

The  residue  series  will  now  be  approximated  by  using  asjmiptotic  expansions  for 


-  S)a  - 


Figure  7 
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fn   and  fg.     The  solutions  of  (l6)  are  asymptotic  to  linear  combinations  of  Hankel 
functions  of  one-third  order.     Noting  that  f,   and  f^  must  satisfy  (l?)  and  (22) 
respectively,  we  find  that 

1/2 


(^3)  ^l-^fe]  4/3(-W) 


ihh)  fg  -  C  f^  ^    j   H0^(_v(x))    ('(n2-a2)l/2i^l](v)  -  iZH^^^Cv) 


% 


Here  the  dot  over  the  Hankel  function  means  differentiation  with  respect  to  its 
argument.  The  functions  v(x)  and  g(x)  are  defined  by 


(U5)     v(x)  =  I  kg3/2  =  k  (    (n2-a2)l/2  ^^ 


=  k  [  (n^. 


From  the  above  definition  we  see  that  g  also  satisfies  the  equation 

(U6)     ggj  =  n^  -  a^. 

In  (U5)j  X  is  a  turning  point  for  (16),  i.e., a  root  of  n(x  )  =  a. 

For  large  k,  the  VJronskian  VJ^Ca)  may  be  obtained  asymptotically  from  (I4.3)  and 
iUh)  •     From  it  we  find  that  (ii2)  becomes 

In  (U?)  we  have  defined  the  operator  jQ.  which  will  prove  to  be  convenient.  From 
the  asymptotic  form  of  V/(a\  t-  can  be  calculated.  Using  this  and  (U3)>  iUh)   in  (lil). 
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we  find 

In  each  teim  of  (U8),   a  is  to  be  evaluated  at  a  ,   a  root  of  (U?). 

We  now  simplify  (U8)  by  utilizing  the  Wronskian  of  the  function  Hi  /*       v(x.  ) 
in  conjunction  xirith  (h7).     We  then  introduce  the  asymptotic  fonn  of  K^/Uv)  and 
H;  /_[y(x  )]   and  carry  out  the  differentiation  with  respect  to  a.     The  expression  (U8) 
for  the  field  then  becomes 

(n^Cx.  )-a^)^/2  expfiJc  J[    +   f  °  |  (n^-a^)      dx  -h  Uca Jy] 

Let  us  denote  the  roots  of  (U7)>   i.e.^the  values  of  v(x.  )  satisfying  (h?),  by 


■- y-     q^'    ,     Now  the  roots  a     are  determined  by  equating  "v(x,  )  to  ■  ,■  ta  e~     q^' 
which  yields 


(50)  kj     (n2-a2)l/2  dx  =  ^  e-i  q^/2  ^ 

We  now  solve  (50)  for  large  k  in  the  same  way  as  we  solved  (19).  The  result  is 


(51)  a^  ^  n(xj^)  +  [6n(x^^)]-^/3  &"^n'(xj^)J^/^  e^"/^  q^t  ... 

From  (5l),  we  can  obtain  the  asymptotic  fonns  of  [nibc,  )  -  a  ]  '  .  If  we  use  this 
and  the  relation  between  the  Hankel  and  Airy  functions  in  (U?)  this  becomes  the 
following  equation  for  a 

(52)  ilV  ^  ^5ni/6  2  r    k     1^/3 


^(%) 


[6n(x,)n.(x,)] 
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The  i\iry  function  in  (52)   is  defined  by 


00 

(53)  A(s)  =  J      cos(t^-st)dt. 


We  next  simplify  the  exponentials  in  (U9).  By  applying  (5l)j  we  find  that 
x^     X 


(5U) 


r-fp- 


2\l/2, 
a  )   '    dx 
m 


a  a 

m  m 


^o       ^^ 


n'-n^^)  -  2.6-1/3  e^"/3[lc-ln(x,)n.(x,)]'\ 


1/2  5) 


dx  +  2 


'.j    (n^-.^)V2a.. 


a 


m 


We  now  substitute  for  the  third  intep:ral  its  value  from  (50)  and  expand  the  integrand 
of  the  first  two  asymptotically  obtaining  four  integrals.     If  the  point  (x,y)  is 
in  the  shadow,  we  can  simplify  the  second  pair  of  integrals  by  means  of  the  equations 
of  the  diffracted  ray  on  which  (x,y)  lies  and  the  limiting  ray  respectively.     This 
yields  the  result 

11/2, 


(55)  {  f .  /  j   (n=-a^)l/2  a.  ~    {f  .  f]  [„=-n^^)j^ 


■dx 


^a      '^a 

m        ra 


^         ^ 


-  [k-ln'(xj^)]^/^  hCx^)]-^/^  e"i/\ly-y(x,)-y(x^^)| 

-  U(3-^  q„)3/2  k-\ 

Here  Y(x,  )  and  y(x,  )  are  respectively  the  ordinates  of  the  points  where  the  limiting 
ray  is  tangent  to  the  boundary  and  where  the  diffracted  ray  splits  off. 

We  now  use  the  relation  between  the  Hankel  and  Airy  functions,  the  differential 
equation  for  the  Aiiy  function  and  (52)  to  rewrite  the  term  ^w,  [v(Xj^)|  /X^AC^^M 
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2/   ^    2 


'     and 


in  (U9).  Inserting  this  result, the  asymptotic  forms  for 

^  v(x.  )  and  the  result  (5^)  into  (U9),  we  obtain  for  any  point  (x,y)  in  the  shadow 

,/^exp(ik{r  +  f°}[n2-n2(x.)ll/WiJcn(x^)lyl 


(,6)  u(x,y)  ^  ^^I^ 


exp 


fe?-''*  q^[6-Wx,)]l/3  k'\)  ly(x^)-Y(x^)|  -l.i(3-\)3/2) 


Here  K(x.  )  is  the  curvature  of  the  incident  ray  at  the  point  of  tangency  to  the 
boundary.  Since  the  second  exponential  has  a  negative  real  part,  the  field  decays 
exponentially  with  the  distance  |y(x.  )  -  Y(x:  )|  from  the  shadow  boundary.  The  result 
(56)  has  been  compared  with  that  of  the  geometric  theory  in  Part  I,  Section  5. 

For  the  problems  corresponding  to  the  limiting  cases  z  ->oo  and  z  =  0,  i.e.,  to 
the  boundary  condition  u  =  0  and  u  =0,  the  eigenvalue  equation  (52)  reduces  to 
^(q  )  =0  and  h   (q^)  =  0  respectively.  In  these  cases  we  can  deduce  the  following 
simple  asymptotic  formulas  for  the  q  for  large  ra. 

2   3/2    ,.36,  ,   j°     "  "° 

In  these  two  cases  we  may  obtain  the  field  u  from  (56)  by  making  use  of  the  equations 
for  q^.  The  reflection  coefficient  R  given  by  (35)  also  simplifies,  in  these  cases, 
to  ^1. 
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U.  The  field  of  a  point  source  in  a  plane  stratified  medium  with  a  plane  boundary 

V7e  now  consider  the  problem  of  a  point  source  situated  at  (x  ,0,0)  in  a  plane 

stratified  medium  with  a  plane  boundary  at  x  =  x,  .  The  index  of  refraction  is  again 

n(x).   The  field  u  is  to  satisfy  (l)  for  x  >  x,  ,  vrith  the  inhomogeneous  term 

5(x-x^)  5(r)  ^ 

o }    (22)  and  the  three  dimensional  radiation  condition.   The  field  will 

tinr 

depend  only  upon  x  and  the  cylindrical  coordinate  r. 

'.fe  define  f(x,a),  the  Hankel  transform  of  u(x,r),  by 


00 

(57)         ^  f(x,a)   =  J    u  jQ(kar)rdr. 


Upon  taking  the  transforr:  of  (l)  we  see  that  f  again  satisfies  (l6)  with  an 
inhomogeneous  term  -  6(x-x  )  as  well  as  (17)  and  (22).   ITnerefore  f  is  given  by 
(26),  as  before.  If  we  take  the  inverse  transform,  we  obtain 

,2  .°° 
(58)    u  =  ^  )  f  Jo(kar)  ada. 

o 
This  can  be  shown  to  be  the  solution  of  the  boundary  value  problem,  provided  that 
f  satisfies  the  above  conditions. 

We  now  use  the  same  procedures  as  in  the  case  of  the  line  source  to  find  the 
asymptotic  form  of  the  field  in  the  lit  and  shadow  regions.  First  the  Bessel 
function  is  replaced  by  its  value  in  teims  of  Hankel  functions.  This  yields  the  two 
integrals 


00 

,  2  , 

(59) 

an 

o 

To  obtain  the  expansion  for  u  in  the  lit  region,  we  split  the  range  of  integration 


in  (59)  into  the  three  intervals  0  to  n(x,  ),  n(x,  )  to  n, ,  and  n,  to  oo.  In  each  of 
these  we  substitute  the  asymptotic  expansions  of  H   '   (kar)  and  the  appropriate 

•t  o 

expansions  for  f  using  (29)  -  (33).   For  the  interval  from  0  to  n(x,  ),  we  obtain  for 
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I^^  an  expression  of  the  form  (3U)  in  which  y  is  replaced  by  r,  the  multiplicative 

•     1  /  ki  V-^^  1/2. 

factor  ri-  by  -  1::?^^  1    and  the  integrand  is  multiplied  by  a   . 
Un    n  1 32nry 

This  integral  can  be  evaluated  asymptotically  by  the  method  of  stationary 

phase,  the  equations  for  the  stationary  points  being  (36)  and  (37)  with  y  replaced 

+  +       + 

by  r.  For  certain  (x,r)  there  are  two  stationary  values  for  I-j^  and  none  for  l^t   I^ 

and  I",  Thus  I*  is  given  by  the  sum  of  the  contributions  from  I^.  If  we  denote 

these  contributions  by  u.  and  u  ,  we  have  u  ~  u ,  +  u  where 


m   u.  ^  ^ 


exp  [ik  I   (n^-a^)"^^^dx  +  ikar] 
(n  -a  )  '^   L"(^o^  "  ®  J 


^o    ^ 


(61)  ^r  ~  H^ 


1/2   ^^l^ll^^/^j 


■  (n^-a2)^/^dx+ikar| 


n  dx 


(XT)^ 


(n2-a2)lA 


[nUf   -  a^j^A 


+     + 

The  remaining  integrals  JZ   and  I"  can  be  evaluated  in  the  manner  outlined  after 

equation  (39) «  For  all  points  in  the  lit  region  the  result  is  the  same  as  that 

just  described  -  u '^  u .  +  u  .  This  result  has  been  compared  with  that  of  the 

geometrical  theory  in  Part  I,  Section  6 . 

Let  us  now  consider  points  x,y  in  the  shadow  region,  /s  before,  we  will 

rewrite  (59)  in  a  way  which  will  turn  out  to  be  convenient  in  this  region.  From 

Figure  7,  we  see  that  an  equation  of  the  form  (Uo)  holds  with  e  ^^  replaced  by 

a  H^  "^(kar).  The  remarks  following  (UO)  concerning  the  case  y  >  0  also  apply  here. 

In  the  same  way  as  before,  we  find  that  u  is  given  asymptotically  by  a  residue 

ika  |y|  Q^ 

series  of  the  form  (UD  vath  e       replaced  by  a  H^  (kar)  and  the  factor  ik 

replaced  by  -i  ik  .  The  a  are  again  the  roots  of  (U2)  which  are  given  asymptotically 

by  (5l).  >fe  utilize  this  result,  (U3),  (lUi),  and  the  asymptotic  expansion  of  the 
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Hankel  functions  to  simplify  the  residue  series.  Then,  using  (52)-(55)  we  further 
simplify  the  result.  We  finally  find  that  for  any  point  (x,r)  in  the  refraction 
shadow  the  field  u  is  given  asymptotically  by 

tx  x^ 
- ,, . ,   .  i^\^\v-^'^^r'^^^^^^^ 

—  ot: —  )  ■? tTTT  I : '5 TTTT — o- PTTT " 


(62)„(x,r)~6-^'^(2|jj 


Here  R(x.  )  and  r(x,  )  are  respectively  the  r  coordinates  of  the  point  of  tangency  of 
the  limiting  ray  (i»e.  the  grazing  incident  ray)  and  the  point  of  shedding  of  the 
diffracted  ray  through  (x,r).  This  result  has  been  compared  with  that  of  the 
geometric  theory  in  Part  I,  Section  6. 


5.  The  field  of  a  line  source  in  a  cylindrically  stratified  medium  with  a  circular 
boundary    ~     — —  —  .  .  _ 

Let  us  consider  a  cylindrically  stratified  medium  for  which  n  =  n(r)  and  which 

is  bounded  internally  by  a  circle  r  =  b.  Let  there  be  a  line  source  parallel  to  the 

z-axis  and  situated  at  (r  ,0).  Again,  the  field  will  be  independent  of  z.  We  require 

that  the  index  n(r)  be  monotonic  and  approach  n,  as  r  tends  to  oo.  For  r  >  b,  the 

field  u  satisfies  (l)  with  an  inhomogeneous  tenn  -r"  6(r-r  )6(&)  on  the  right  side. 

On  the  boundary  u  must  satisfy  (22)  with  ~  replaced  by  ~  and  at  infinity  the 

dx  or 

radiation  condition  (23) » 

Since  u  has  to  be  even  and  periodic  in  (^,  we  assume  that  the  field  can  be 
written  in  the  foim 
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(63)      ^  =  ■^  Z   ^^^''"^ 


2n 

rri=-oo 


e   . 


The  function  f(r,m)  is  expressible  as  the  Fourier  integral 

2n 
(6U)       f(r,m)  =  J   ue~    d9. 


We  now  apply  this  finite  Fburier  transform  to  (l)  with  the  right  side  properly- 
modified.  We  find  that  the  function  r    f(r,in)  satisfies  (16)  with  a  replaced 

O   T        ^O  — 1/2 

by  (m  -  T-)(kr)''"  and  with  the  inhomogeneous  tenn  -r    5(r-r  )  on  the  right  side. 

The  solution  to  this  equation  satisfying  (l?)  and  modified  (22)  is  given 
by  an  expression  of  the  form  (26).  However,  in  it  x^,  x^,  and  a  must  be  replaced 
respectively  by  r^^,  r  and  m  and  the  entire  expression  must  be  multiplied  by  r~  . 
With  this  value  for  f(r,m),  we  can  verify  that  (63)  is  the  solution  to  our 
boundary  problem.  However,  to  obtain  an  asymptotic  expansion  for  u,  it  is 
convenient  to  rewrite  (63)  by  using  a  Watson  transformation. 

By  Cauchy's  residue  theorem,  we  can  transform  the  sum  in  (63)  into  the 
contour  integral 


(65)      u=        1    J        t^:^^^ da. 


C,^C2 


Here,  CL  and  C^  are  the  two  straight  lines  Im  a  =  ±  constant  in  the  neighborhood  of  the 
real  axis.  The  direction  of  integration  along  C^  is  from  left  to  right  and  for  C^  from 
right  to  left.     Using  the  fact  that  f  is  even  in  c  we  can  write  (65)  as 

^    f(r,<i)cos  a(Q-n) 
(66)         ---■h)  sinna  ^°- 

^1 
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To  put  (66)  in  a  form  which  v;ill  yield  the  field  in  the  lit  region,  we  set 
cos  a(&-n)  =-ie   sin  na  +  e   cos  a©.  We  then  move  the  contour  for  the  first 
summand,  which  now  has  no  poles  on  the  real  axis,  dovm  to  the  real  axis  and  then 
let  a  =  ka  in  both  integrals.  As  a  consequence,  we  have 

(67)   u  =  ^  rf(r,ka)ei^da  -  ^  f  f(r,ka)e^"^-cos  akQ  ^^  ^ 
2n  ^     '  2tti  j  sxn  Tiak 

-00  C- 

The  first  integral,  being  of  the  form  (27)j  can  be  treated  as  in  the  previous  two 
sections.  If  we  divide  the  positive  a-axis  into  the  three  intervals  0  to  n(b). 


n(b)  to  a,  and  n^  tooo,  we  find  that 


dr+ika&l 


n(b)   s'^pl      P-  T        dr+ik^+R  expjik  1    f  +  |    >(n^-  %  | 

(68)    i^-^r       ^-^^'^    "^  /       \   i^   Ml    ^i y^^ 

^  ^^  J  r  2  2     211/1,  r  ?  ?,   .     ?ii/). 


Here 
(69) 


[r^n^  -  a^^A  [r^n^(r,)  -  a^]^A 

r  2  2/,  N   2I1/2  ^ ,  ^ 

Lb  n  (b)-a   J         +  bZ 


\bV(b)-a^" 


^-bZ 


We  now  apply  the  raethod  of  stationary  phase  to  I,  and  find  two  stationary  points 
for  certain  values  of  (r,9).  The  integrals  L,  p  and  I~,  for  these  values  of  (r,©), 
do  not  have  stationary  points.  The  contributions  to  the  field  by  the  first  integral 
of  (67)  is  thus  the  sura  of  the  contributions  from  the  two  stationary  points  of  L., 
If  we  designate  these  contributions  by  u.  and  u  ,  we  get  the  geometric  field 


u^  -%/  u.  +  u  where 
g     1    r 


exp  I  ik  I  (n^  -  2_)   dr  +  ik^  +  ^" 


^2  1/2 
(70)   Uj_ ^ ^-- ^—- -i72- 


\  2  2f      X  2"ll/Uf  2  2  2]l/h(a,     f  ^  rTLdr__ 
LV  ^^o)-^  J  ^  b  n  -a  J  ^  (^8nk  J  ^-^-^--^-^ 
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V2 


(71)     u„   ^ 


K-i'le-^j 


R   exp  I  ik  1  I   +   1      (*  (n"  -  ~  J      dr  +  ika9  +  ^" 


r  r  '  2_/2 

The  remaining  integrals  !<,  -,,   I~  are  handled  as  in  Section  3  and  similar  results 
are  obtained. 

Vie  now  show  that  the  second  integral  of  (67)  :\T.elds  the  diffracted  field  u,  . 
Vfe  first  make  the  change  of  variables  a  =  ka.     Then  we  take  the  resulting  integral 
over  a  finite  segment  of  C   .     This   is  equivalent  to  an  integral  along  a  semicircle 
in  the  upper  half  plane  plus  the  residues  of  the  integrand.      Again,  it  can  be  shown 
that  the  integral  over  the  semicircle  tends  to  zero  as  its  radius  approaches 
infinity.     We  are  left  with  the  residue     series 

e^"^^m  cos  ka  e     f, (r^,ak)f5(r^,ak) 
(72)     u^^^iL':^  ml>  d     < 


d       ^o     m  |-W(a   ) 

8a       m 


Here,  the  a  are  the  zeros  of  the  Wronskian  W(a). 
m 

The  analysis  now  proceeds  along  previous  lines.  The  eigenf unctions  f-,*i"p 
are  replaced  by  their  asymptotic  representations  (U3)j  (UU).  However,  the  function  v 
is  now  defined  by 


r 


^a 


2  1/2 


(73)    v(r)  =  I  (n^  -  ^  )    dr. 


2  2  2 

The  quantity  r  is  a  zero  of  r  n  -a  and  corresponds  to  a  turning  foint.  Instead 

2    2  2 
of  (U6),  g  now  satisfies  the  equation  gg  =  n  -a  .  Then,  asymptotically  for  large  k 

r  p      p  j  ]  /p         12    ?  "2 1 1  /2 

the  zeros  a  are  given  by  (hi)   with  In  (x,  )  -  a  j  '  replaced  by  In  (b)-a  b  J  ' 
and  with  v  defined  by  (73).  These  eigenvalues  can  be  approximated  as  in  Section  3. 
Instead  of  (5l),  we  obtain 

C7U)   a^  ^.n(b)  .  V-/3  [^f/f-r^lS^^Y/^  . 
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The  q  satisfy  (52)  mth  the  factor  n"   (x^)  jn' (x^)J~    replaced  by 

We  no'^  can  carry  out  approximations  and  simplifications  similar  to  those 
foUovdng  equations  (U?)  and  ($l)  in  Section  3.  As  a  consequence  we  obtain  for  the 
diffracted  field  at  any  point  (r,9),  the  result 


(75)       u,(r,&) 


ne 


in 


/3  r 


U-6 


Tn 


b  K(b)-»-l 
kbn(b) 


2/3 


exp 


')m;jF'-^ 


(b) 


Lb  n  (b)         J  LbV(b) 


1/2 


dr 


TTT 


m 


exp[-Ui(3"\)^^^] 

[l-exp[2Tiikbn(b)+2ne5"^/^S^]j'qJ.^(qj^)+3A'^(q^)] 


•  I   exp[ikbfl(b)(2Ti40)  -^  e^"^/^q^(6-^bn(b)}^/3[bK(b)+lJ^/5(2Ti  +  ©^  -  Y^ )] 


+  exp 


iicbn(bX2n-  e)  +  e^"^/\{6~^bn(b))^/^[b,K(b)+4^/^(2n-2e  +  9^-  I^^ 


I      In  this  expression,     K  (b)  is  the  curvature  of  the  limiting  ray  at  the  point  of 

tangency  (b,f,  )  to  tlie  circular  boundary.     (b,e^)  is  the  point  on  the  boundary  where 
the  diffracted  ray  splits  off.     The  quantity  S     is  defined  by 


(76)       S^  =  q^[6-^bn(b)]^/3   [bK(b)  +  lj2/^. 


Another  fomula  for  the  diffracted  field  can  be  derived  from  (66)  in  exactly 
the  same  way.  We  find  instead  of  (75) 
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r   r. 
ni/3  ^KW.^..  12/3  ^^'^ 


TTT 


m 


+  exp&bn(b)(2n-  &)  +  e^"""^^  q^(6"W(b)f  ^/^/b,K  (b)+l}  ^/^(2ti-294^^-¥^) 

This  form  is  valid  in  the  refraction  shadow  whereas  (75)  holds  in  the  lit  region. 
The  two  sums  in  each  expression  are  the  modes  of  the  cylinder.  From  (77)j  we  can 
get  the  result  (56)  of  Section  3,  /ie  set  r  =  x+b,  r  =  x  +b,  9  =  b"  y  and  let  b 
tend  to  00.  On  the  other  hand,  by  settinr,  n'  =  0,  we  obtain  the  asymptotic  form 
of  the  diffracted  field  due  to  a  line  source  in  a  homogeneous  medium. 

The  results  (70)j  (71)  and  (77)  have  been  compared  with  the  geometrically 
constructed  fields  of  Part  I, 
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APPENDIX 

The  mode  expansion  of  the  field  due  to  a  line  source  in  a  plane 
stratified  medium  with  a  plane  boundary  is  obtained  by  applying  a  Watson 
transformation.  We  now  show  that  the  branch-cut  integral  and  integral  along 
the  pure  imaginary  axis  which  arise  in  connection  with  the  transformation  are 
asymptotically  zero.  Also,  the  integral  along  the  semi-circle  tends  to  zero 
with  increasing  radius. 

By  using  (U3)  and  (UU)  we  can  express  the  integrand  of  (28)  asymptotically 
for  large  k.  We  then  obtain  for  the  integral  the  expression 


(Al) 


A  similar  result  involving  e"  ^^  holds.  However,  the  branch  cut  on  the  real 
axis  occurs  for  only  one  of  these  expressions,  depending  on  how  we  close  the 
contour  in  the  complex  a-plane. 

Due  to  the  raulti-valuedness  of  v(x  )  (or  v(x)  and  ▼(^c^))  its  argument 
will  differ  on  the  opposite  sides  of  the  cut  B  .  Consider  first  v(x^).  We 
subdivide  the  real  positive  axis  into  the  intervals  0  to  n(0),  n(0)  to  nCx^), 
n(x  )  to  n.   and  n,  to  oo.  Then  values  correspond  to  particular  turning  points 
of  (16),  Part  n.  Now  the  value  n(x  )  is  a  branch  point  of  v(Xq).  Because  of 
the  monotonicity  of  n(x)  (see  Figure  la),  it  follows  thai  to  every  value  of  a 


*The  cut  actually  extends  fix)m  -n(x^)  to  n(x^)  on  the  real  axis.  For  the 
analysis  of  (Al),  we  consider  only  the  part  on  the  positive  real  axis. 
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in  the  interval  n(0)  to  n,  corresponds  a  real  turning  point  x^.  In  particular 
for  a  in  the  interval  n(0)  to  n(x^),  the  turning  point  x„  <  x^.  Therefore, 

the  quantity  (n  -  a  )  '  is  real  and  the  argument  of  vCx  )  is  zero.  For  those 

2   2  1/2 
a  in  the  interval  n(x  )  to  a,,  the  turning  point  x^  >  x^  so  that  (n  -  a  ) 

2   '  3/2 

is  pure  imaginary.  Near  the  turning  point  v  ~  »  n  (x^)(x  -  x^^)  '  .  Thus, 

the  argument  of  v(x  )  for  these  values  of  a  is  -3't/2. 

Corresponding  to  the  intervals  (0,n(0))  and  (n-,00)  there  are  no  real 
turning  points  so  that  the  argument  of  v(x  )  in  the  first  interval  is  zero  and 
in  the  second  -3n/2.  On  the  other  side  of  the  cut,  its  argument  frcm  0  to  n(x^) 
must  be  -3n. 

In  the  same  way,  we  can  show  that  the  arguments  of  v(x)  and  vCx,  ) 
are  zero  and  -3n  on  either  side  of  the  branch  cuts  determined  by  their  respective 
branch  points  of  n(x)  and  n(x,  )  and  -3«/2  for  values  of  a  beyond  tho  branch 
point. 

With  the  previous  remarks,  we  find  that  the  expression 

has  on  one  side  of  the  branch  cut  up  to  the  branch  point  n(x^)  the  argument 
zero  and  on  the  other  side  the  argument  -2n.  For  values  of  a  >  n(x^),  its 
argument  becomes  -  n/2  and  -3^/2  respectively. 
By  making  use  of  the  relations 


we  now  can  show  that  (Al)  is  zero.  In  the  interval  (n(x,  ),  n(x  )),  for 
example,  the  two  products  resulting  from  the  multiplication  of  (A2)  with  the 
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Hankel  fimctions  of  argrments  v(x)  and  v(x  )  aiB  positive  on  one  side  of  the 
cut.     Since  on  the  other  side,    (A2)  has  an  argtunent  of  -2n  and  then  Hankel 
functions  have  arguments  of  -3n,  application  of  (/3)   and  {Ah)  shows  that  these 
products  are  still  positive.     On  the  other  hand  the  argument  of  _[\'^ /i  /_Q,^\/i 
doesn't  change  in  this  interval.     Since  the  paths  of  integration  are  traversed 
in  opposite  directions,  the  integral  over  this  interval  is  zero.     In  the  same 
way,  we  can  show  that  (/,l)  vanishes  over  the  remaining    portions  of  B. 

On  the  imaginary  axis  the  argument  of  v  changes  from  zero  to  -3n  in 
going  from  Ima>0tolma<0  and  that  of  (A2)   changes  from  zero  to  -2n, 
Therefore,  if  we  set  a  =  -iz  and  utilize   (A3)  and(AU),  we  find 


(AS) 


OQ 


OQ 


W^W¥(Z>^ 


H()^[v(x,)e-3"i]U)[v(^,e-3"q  .  y^y^l.^^U'^']^ 


(2) 
■3 


n.^% 


r(x,)e-3"i[ 


^^V^ 


dz 


k 


-  4)^^=^'!  ^f5^ 


^     -  f    f  e^^y  dy. 


dz 


Thus,  the  sum  of  the  two  integrals  along  the  imaginary  axis  is  asymptotically 
equal  to  zero. 
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If  we   nov;  expand  the  Hankel  functions  in  (Al)  asymptotically  for  large  \a\, 
we  find  that  the  integrand  behaves  like 

-k(3^-x^)p  cos  9  -k(2x^-x^-x)p  cos  & 

(a6)    e  +  Re 

In  (a6),  a  =  pe^®  and  R  =  0(1)  for  large  p.  Along  17  j  (^  S  arg  a  <  J  so  that 


(A7) 


But 


/   ^1  /'^^  -k(x.-x^)p  cos  &       -k(2x, -X  -x)p  cos  e\ 
J   '0{  )       V  d9  +  Re     "  °         )d& 


//^  -k(x^-x^)p  cos  9     ^  r    -kp(x -X  )1 
(A8)    je      ^        d^<^l-e        J- 

o 

For  X  /  X  ,  this  integral  tends  to  zero  as  p  ->oo.  In  the  same  way,  the  second 
integral  of  (A?)  can  be  shown  to  approach  zero.  When  x  =  x^,  the  term  corres- 
ponding to  the  first  exponential  of  (A?)  reduces  to  unity.  However,  instead  of 
using  the  bound  of  unity  on  |e  ^^| ,  we  estimate  its  integral  and  this  yields 
an  expression  given  by  an  equation  such  as  (A8),  As  a  result,  the  integral 
along  I^  tends  to  zero  for  all  x  as  p  ->oo. 

Finally,  we  note  that  the  analogous  manipulations  for  the  problems  of 
Sections  h   and  5  of  Part  II  may  be  justified  in  the  same  way. 
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